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ABSTRACT 


A  theoretical  study  is  made  of  hypersonic  flows  over  thin  bodies  subject 
to  the  combined  effect  of  leading-edge  bluntness,  boundary-layer  displacement  as 
well  as  surface  incidence .  Within  the  scope  of  the  boundary-layer  theory  and 
the  hypersonic  small-disturbance  theory,  an  approach  is  developed  based  on  a 
flow  model  consisting  of  three  adjoining  regions?  an  inner  laminar  boundary 
layerj  an  outer  (detached)  shock  layer j  and  between  these  two,  a  low-density 
entropy  layer »  This  approach  allows  simultaneous  treatment  of  the  boundary  layer 
and  its  outer  invlscid  flow  field  in  a  manner  consistent  with  the  Newtonian 
shock-layer  approximation,, 

The  theory  pertaining  to  the  leading  approximation  for  strong  shock  and 
y-^i  is  completely  developed  and  its  simplicity  permits  solutions  to  a  number 
of  problems  of  interest.  For  the  flat  plate  at  zero  incidence,  the  theory 
yields  a  solution  agreeing  with  the  blast-wave  theory  at  one  limit  and  the 
strong  inviscid -viscous  Interaction  theory  at  the  other.  Results  are  also  ob¬ 
tained  for  the  flat  plate  at  incidence.  The  results  for  the  pure  inviscid  flows 
with  tip-bluntness  effect  exhibit  an  oscillatory  decay  in  pressure  on  slender 
cone  and  wedge  afterbodies.  This  phenomenon  signifies  an  idealized  form  of  re¬ 
impingement  of  the  shock  layer  on  the  afterbody  at  glancing  incidence. 

Complementary  to  the  aforementioned  study,  an  examination  is  made  of  the 
similitude  of  hypersonic  boundary  layers,  taking  specifically  the  leading-edge 
bluntness  into  account.  Rules  for  correlating  pressure,  shock  shape,  skin  fric¬ 
tion  and  heat-transfer  rate  are  given.  The  simplification  of  these  rules  for 
flows  with  strong  shock  waves  is  pointed  out. 
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pressure 

Prandtl  number  ^  ^ 
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Reynolds  number 
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temperature 
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respectively 

2  ^  the  variables  related  to  ^  and  /p  ,  respectively,  as  defined 
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the  angle  of  attack  of  the  flat  plate »  Positive  oC  refers  to 
compression  or  wedge  flow  over  the  surface  in  question^  also  half¬ 
cone  angle 

,  specific-heat  ratio 

the  boundary-layer  thickness  and  displacement  thickness,  respec¬ 
tively 


respectively 


the  shock  angle,  i»e,  , 

flow-deflection  angle  ^ 


which  may  be  taken  as  the  typical 
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Q 

K 

A 

I)r) 

TT 

r 


cd 


H-H^ 

a  function  whose  derivative  d§/9fl 
thermal  conductivity  of  the  gas 
viscosity  of  the  gas 
certain  variables  in  Section  1| 

two  spatial  variables  used  in  the  boundary-layer  analysis,  see 
Eq«  (2ol) 


a  number  which  is  zero  for  plane  flows  and  unity  for  axisymmetric 
flows 


gas  density 

Ke/^; 

a  body  thickness  parameter  used  in  the  similitude  consideration 
(for  the  wedge  and  cone,  s  ) 

€  (0,661;  +  1<,73  -5^^)  y  f  a  parameter  governing  the  boundary- 
layer  displacement  effect 

the  constant  exponent  in  the  relation  jUCQ  f 


Subscripts 
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e 

L 
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pertaining  to  the  inner  edge  of  the  entropy  layer 
pertaining  to  the  Blasius  solution 
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based  on  the  reference  length  ~  L') 

pertaining  to  the  condition  behind  the  strongest  part  of  the  shock 
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NOTATION  (Contd) 


behind  the  shock  (at  /jC  ) 

based  on  the  nose  thickness  (/^st) 

at  the  body  surface  or  wall 

at  the  outer  edge  of  the  boundary  layer 

pertaining  to  the  free-stream  stagnation  conditions  (for  the  tem¬ 
perature) 

pertaining  to  the  free-stream  condition 
based  on  the  reference  temperature 


LIST  OF  ILLUSTRATIONS 


Figure  Page 

1(a)  Illustrations  of  the  Streamline  Pattern  of  the  Flow  Field 

Around  a  Slender  Body  mth  Small  Blunt  Nose  .........  6? 

1(b)  Division  of  the  Inviscid-Flow  Region  About  a  Slender 

Afterbody  Into  a  Shock  Layer  and  an  Entropy  Layer  o  .  .  «  .  6? 

2  Determination  of  the  Maximum  Perturbation  Velocity 

|Au|  3=  ui.  Over  a  Flat-Plate  Afterbody  in  an  Inviscid 

Hypersonic  Flow  ,  o  o  .  o  «  .  .  .  .  o  .  »  »  «  »  .  «  «  «  «  »  68 

3  Results  Based  on  the  Zero-Order  Theory  of  a  Flat  Plate 
Under  the  Combined  Effect  of  Tip  Bluntness  and  Boundary- 

Layer  Displacements  o  ......  o  o<>  .<.0  «  69 

ii  Results  Based  on  the  Zero-Order  Theory  of  a  Flat  Plate 

Under  the  Combined  Effect  of  Boundary- Layer  Displacement 

and  Angle  of  Attack  o9«ooooooooooooooooo  TU 

5  Results  Based  on  the  Zero-Order  Theory  of  a  Flat  Plate 
Under  the  Combined  Effect  of  Leading-Edge  Bluntness  and 

Angle  of  Attack  ooooo99»ooooeoo  9999999  T1 

6  Results  Based  on  the  Zero-Order  Theory  of  a  Blunted 
Slender  Cone  909999999999999999999999 

7  Results  Based  on  the  Zero-Order  Theory  of  Hypersonic  Flow 

Over  Thin  Wedge  Subject  to  Leading -Edge  Bluntness  and 
Boundary-Layer  Displacement  Effects  «  9  9  ,  .  9  .  «  9  .  .  .  73 


X 


SECTION  1 


INTRODUCTION 


Two  outstanding  features  characterize  the  hypersonic  flow  around  a  thin 
body*  One  is  the  displacement  effect  of  the  boundary  layer  and  the  other  is  the 
downstream  influence  of  a  small,  but  blunt,  nose  or  leading  edge*  The  pro¬ 
nounced  displacement  effect  is  related  to  the  viscous  dissipation  of  the  high¬ 
speed  flow,  which  raises  the  local  temperature  and  thus  the  kinematic  viscosity. 
The  lateral  extent  of  the  hypersonic  boundary  layer  may  therefore  become  com¬ 
parable,  or  greater  than,  the  thickness  of  the  thin  body,  and  an  interaction  of 
the  inviscid  and  viscous  flow  regions  must  take  place .  The  importance  of  the 
leading-edge  bluntness  may  be  perceived  from  the  connection  between  the  energy 
gained  by  the  cross-flow  field  and  the  nose  drag.  Since  the  drag  of  a  blunt 
nose  is  rather  high,  the  work  done  by  it  alone  may  account  for  a  sizable  portion 
of  the  total  energy  delivered  to  the  fluid.  The  unsteady  motion  observed  in  the 
cross-flow  plane,  that  is,  a  plane  fixed  with  respect  to  the  ambient  fluid  and 
transverse  to  the  direction  of  motion,  may  therefore  be  compared  with  that  of  a 
blast  wave  generated  by  instantaneous  release  of  energy  of  an  amoimt  proportional 
to  the  nose  drag.  For  hypersonic  flight  of  slender  bodies,  leading-edge  blunt¬ 
ness  is  generally  essential  for  heat-transfer  control,  while  for  flight  at 
relatively  high  altitude,  the  boundary-layer  displacement  effect  is  important. 
This  study  will  consider  hypersonic  flows  in  which  both  of  these  phenomena  are 
present. 

There  is  a  body  of  theory  dealing  with  each  of  these  two  effects  as  inde¬ 
pendent  phenomena.  For  the  boundary-layer  displacement  problem  taken  separately, 
we  have  the  rather  extensive  works  on  strong  and  weak  interactions  of  the  invis- 
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cid  and  boundary-layBr  flows  over  a  (shai^j-edged)  flat  plate,  as  well  as  treat- 

1-17  ' 

ments  of  related  problems  »  These  theories  provide  an  account  of  the  self- 
induced  pressure  in  the  outer  flow  associated  with  the  displacement  effect  of 
the  boundary  layer.  From  the  literature  on  the  subject,  the  early  contributions 
of  Shen^^^,  Lees  and  Probstein^®^,  Li  and  Nagamatsu^,  and  also  the  definitive 

9 

work  of  Stewartson  should  be  noted.  Taking  the  bluntness  effect  separately,  we 

l8  19  20-2li 

have  the  strong  blast-wave  theories  of  Lin  ,  Lees  and  Kubota  and  others 

These  works  make  use  of  the  analogy  between  the  transient  cross-flow  field  down¬ 
stream  of  a  blunt  nose  and  the  flow  field  produced  by  a  violent  explosion  of  the 

2^  26 

type  previously  treated  by  Taylor 

On  the  other  hand,  for  the  problems  in  which  both  the  displacement  and 

bluntness  effects  are  important,  there  has  not  been  any  satisfactory  theory.  A 

theory  accounting  simultaneously  for  both  effects  is  essential  not  only  for  its 

practical  Implication,  but  also  because  the  task  of  isolating  these  two  effects 

27-31 

in  hypersonic  flow  is  generally  not  simple,  even  under  laboratory  conditions 
The  aim  of  this  paper  is  to  present  a  simplified  but  consistent  approach  to  the 
problem  which  will  render  possible  the  study  of  interaction  of  the  boundary  layer 
with  the  outer  flow  under  the  influence  of  leading-edge  bluntness  as  well  as 
afterbody  geometry. 

A  continuum  flow  model  consisting  of  a  laminar  boundary  layer  adjacent  to 
the  slender  body  and  an  inviscid  outer  flow  field  will  be  assumed  for  this  theo¬ 
retical  study.  Since  the  Mach  numbers  considered  are  high,  and  the  bodies  thin, 
the  lateral  extent  of  the  flow  field  bounded  by  the  bow  shock  is  small.  Our 
analysis  may  therefore  be  carried  out  within  the  scope  of  the  boundary- layer 
theory^^^  and  the  (inviscid)  hypersonic  small-disturbance  theory^^”^^.  Within 
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this  framewrlc^  tm>  complementary  studies  are  made*  Ilrstj  and  the  major  part 
of  our  analysis,  is  the  development  of  an  approach  -ssihich  takes  account  of  the 
combined  effect  of  boundary- layer  displacement,  leading-edge  bluntness,  and  body 
thickness.  To  achieve  this  end,  certain  approximations  and  simplifications  ■sdLll 
be  introduced.  The  second  study  is  concerned  -with  hypersonic  similitude  •ciien 
both  boundaryi-layer  and  tip-bluntness  effects  are  present,  and  mth  simplifica¬ 
tion  of  the  similarity  law  for  flows  with  strong  shock  waves.  For  this  simili¬ 
tude  consideration,  a  number  of  the  assumptions  introduced  in  the  first  study  are 
not  required. 

Underlying  the  present  ^proach  are  the  two  basic  concepts;  a  local  simi¬ 
larity  in  the  boundary  layer  and  a  detached  shock  layer  in  the  outer  inviscid 

flow.  Both  of  these  concepts  are  desirable  for  a  development  consistent  with  the 

3^-il.O 

so-called  Newtonian  shock-layer  approximation  ’  #  The  consistency  between 

local  similarity  and  the  shock-layer  approximation  has,  in  fact,  been  discussed 
ill 

earlier  ,  The  concept  of  a  detached  shook  l^er,  i^ch  includes  the  so-called 

39, ill 

"free  layers”  as  special  cases,  is  a  necessary  modification  of  the  shock- 

layer  theory  to  allow  treatment  of  tip  bluntness, 

A  similar  treatment  of  the  tip-bluntness  effect  has  been  given  previously 
21  22 

ly  Chernyi  *  ,  The  essential  idea  of  a  detached  shock  layer  is  Implicit  in  his 

integral  method  for  iypersonic  inviscid  flows,  "While  Chernyi 's  method  and  the 
present  analysis  are  equivalent  in  the  idealized  limit  •/ — ,  they  differ 
in  the  higher  order  approximation  for  Y  not  close  to  unity*  It  m^  also  be 
pointed  out  that  the  inviscid  tip-bluntness  effect  has  been  studied  recently  by 
Guiraud^^^^^,  Here,  the  idea  of  a  detached  shock  layer  was  eniployed,  and  certain 
equations  obtained  are  very  similar  to  those  to  be  developed  later  in  this  re- 
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port.  Guiraud^s  study  does  not  lead  to  an  explicit  solution  of  the  problaa^ 
howeTera 

Throughout  the  present  analysis  an  ideal  gas  mth  constant  specific  heats 
id.ll  he  asstoaed  ahead  of  ajid  behind  the  shock  (althou^  y  may  change  across 
the  shock)  o  In  order  to  arrive  at  a  simple  form  of  local  similarity  in  the 
boundary  layer the  Ih’andtl  number  is  taken  to  be  unity,  the  bo<^  temperature 
is  assumed  to  be  uniform,  aid  a  linear  vd.SGosity->temperature  relation  is  adopted. 
The  present  study  is  primarily  concerned  -with  plane,  steady  floira.  For  the 
purpose  of  comparing  plane  and  axisymmetric  flews  for  the  inviscid  tip“bluntness 
effect,  however,  the  problem  of  a  blunted  slender  cone  is  also  studied  along 
mth  the  corresponding  problem  of  plane  flow  over  a  blunted  wedge. 

The  three  requirements  which  control  the  accuracy  of  the  followi.ng  analysis 
of  the  present  flow  model  may  be  expressed  as  follows? 


(a)  The  disturbance  velocity  ss  \~ —  /  ,  as  well  as  the 

U  \JJ  J  ' 

square  of  the  flow  angle,  in  the  inviscid  flow  field  must  be 
negligibly  small  in  comparison  to  unity. 

(b)  The  bow  shock  wave  must  be  sufficiently  strong. 

(c)  The  specific  heat  ratio  /  must  be  sufficiently  close  to  unity* 

The  first  requirement  implies  eiclusion  of  the  blunt-nose  region  itself 

fitsm  the  present  enquiry,  and  we  will  deal  only  with  its  downstream  effect.  The 
consistency  of  these  requirements  and  the  degree  of  the  approximation  involved 
will  be  examined  later.  The  aero-order  theory,  i&ich  is  to  be  studied  in  detail. 


The  axisymmetric  problan  involving  boundary-llayer  displacement  is  complicated 
by  the  transverse- curvature  effect.  In  principle,  however,  a  modified  form  of 
the  local  similarity  can  be  obtained  for  the  boundary  Isyer,  although  it  ii\ay 
lose  the  appeal  of  simplicity. 
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pertains  to  the  leading  approximation  under  Ihe  last  tw  requirements  cited^  and 

35»i|.0 

may  therefore  be  compared  mth  the  classical  result  of  BuseiE^m  (NewboiuLan- 
Plus-Cent3rifugal“Force)  o 

In  -view  of  the  continuum  assumption j  the  present  theoretical  model  cannot 

be  used  as  a  valid  description  of  flow  phenomena  idien  the  molecular  nature  of 

the  gas  is  manifestedo  On  the  other  hand^  the  effects  of  velocity°slip,  and 

temperature- Jump,  as  well  as  of  the  "vorticity  interaction”  -sdiich  precedes  the 

breakdown  of  the  usual  continuum  model,  may  be  seen  to  belong  to  the  same,  or 

13,l6,ijlj. 

higher,  order  as  the  error  in  the  standard  boundary-layer  approximation  e 

¥e  may  th^efore  expect  the  continuum  theory  to  remain  meaningful  so  long  as 
there  exists  a  sufficiently  thin  boundary  layer,  and  this  condition  is  implied 
by  requlremmit  (a)*  Based  on  the  solutions  obtained,  it  is  possible  to  define 
a  domain  of  validity  for  the  present  continuum  flow  model  in  terms  of  the  Mach 
and  Reynolds  numbers  of  the  free  stream*  It  may  be  seen  that  this  domain  is 
large  enough  to  permit  many  useful  applications*  The  adoption  of  the  unit 
Prandtl  number  and  a  linear  viscosity-temperature  relation  may  seem  at  first 
sight  to  be  oversimplifying  the  problem*  As  examples  will  reveal,  these  sim¬ 
plifications  are  not  likely  to  affect  the  skin  friction,  heat  transfer,  and  the 
displacement  thickness  in  a  significant  way,  particularly,  •tftien  an  appropriate 
reference  temperature  is  chosen  for  the  viscosity-temperature  r^ation*  ¥e 
may  notice  also  that  for  the  present  analysis,  these  two  assumptions  are  not 
strictly  necessary*  Analyses  similar  to  the  present  one  may  be  carried  out  for 
an  arbitrary  but  constant  Prandtl  number  and  a  more  general  viscosity  law  with 
no  apparent  difficulty*  Another  limitation  of  the  present  study  is,  of  course, 
the  assumption  of  constant  specific  heats*  The  gas  behind  the  normal  shock  in 


front  of  a  Mtmb  nose^  under  most  circumstances^  mil  be  in  a  dissociated^,  or 
internally  excited  states  the  specific  heat  ratio  of  the  gas  isill  be  closer  to 
unity  than  that  under  standard  atmospheric  conditiono  Generally  speaking,  the 
shock-layer  approximation  may  Improve  as  a  result  of  the  real-gas  properties  at 
high  terapecatureo  Stilctly  speaking,  however,  the  specific  heat  ratio  or  its 
equivalent  cannot  generally  be  taken  as  constant  for  the  entire  flow  fieldo  If 
the  flow  field  over  the  slender  afterbody  can  be  assumed  to  be  near  local 
theimal  equilibrixim,  and  if  the  body  temperature  is  sufficiently  low,  the  real- 
gas  effects  may  not  be  critical  since  the  local  temperature  downstream  is  con- 
siderably  lower  than  the  stagnation  temperature  at  the  nose«  In  any  case,  the 
approach  presented  herein  may  still  be  enqiloyed  to  simplify  the  analysis,  al¬ 
though  it  may  not  lead  directly  to  a  solution*  This  point  will  be  discussed 
further* 

A  somewhat  shortened  version  of  this  investigation,  together  with  exparl- 
mental  results,  has  been  presented  in  Eef*  i|5  *  In  addition  to  the  inclusion  of 
tip-bluntness  effects  in  axisymmetric  flows,  further  discussion  and  clarifica¬ 
tion  not  given  in  Ref*  will  be  taken  up  here* 

In  the  following  sections,  the  local  similarity  in  the  boundary  layer  and 
its  consistency  with  the  inviscid  shock-layer  approximation  -wLU  first  be  dis¬ 
cussed*  The  theoiy  of  a  detached  shock  layer  will  next  be  considered*  The 
zero-order  theory  for  the  complete  problem  will  then  be  developed  and  applied 
to  particular  cases  of  interest* 

^  The  experimental  study  described  in  Ref*  ii5  153,3  conducted  in  the  Cornell  Aero¬ 
nautical  Laboratory  11  x  1^-inch  Hypersonic  Shock  Tunnel  under  the  sponsorship 
of  the  TJ*  So  Air  Force  Office  of  Scientific  Research* 
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SECTION  2 


THE  LOCAL  SBIELARITr  AND  THE  SHOCK-»LAIM  APPROnMATION 

The  idea  of  assuming  a  local  similarity  in  the  boundary  lasrer  as  a  means 
of  estimating  heat-» transfer  rates  on  blunt  bodies  at  high  speed  -was  first  proposed 
by  Lees^^  and  utilized  by  others^^  Lees  observed  that  the  enthalpy  distri¬ 
bution  in  the  boundary  layer  is  rather  insensitive  to  pressure  gradient 5  parti¬ 
cularly  if  the  boundary  layer  is  cooledo  For  most  engineering  appli cation j,  it 
suffices  to  ignore  the  pressure  gradient^  except  for  its  effects  on  changing  the 
local  Reynolds  numbero  In  this  case  the  problem  is  reducible  to  that  of  a  flat 
plate^,  and  for  this  problem^  a  self-similar  solution  exists  if  the  "wall-to- 
stagnation  temperature  ratio  is  constant  along  the  plat©o 

It  is  essential  to  recognize  that  in  the  present  study 5  ^ere  only  thin 
bodies  are  considered^  the  pressure  gradient  is  unimportant  not  primarily  because 
the  wall  is  made  cold^  but  because  the  direct  effect  of  pressure  gradient  on 
the  hypersonic  boundary  layer  is  at  most  of  the  order 

y-i 


The  condition  of  anall  e  is  a  requirement  in  listing  shock-layer  theories  for 
inviscid  hypersonic  flows  o  It  is  riatural  to  treat  the  present  boundary- 
layer  problem  by  an  analogous  perturbation  procedure  which  uses  e  as  a  control¬ 
ling  parameter*  The  argument  used  for  the  local  similarity  method  in  blunt  body 
applications  is  avoided  in  the  present  approach* 

'' 

In  a  more  specific  sense^  the  similarity  under  discussion  should  be  termed 
the  "local  flat-plate  similarity"* 


? 


The  foregoii^  considerations  are  readily  seen  after  transforming  the  var¬ 


iables  /)[/  s  ^  xn  the  boundary  layer  as  follows  s 

s-f':  i 


L 


n 


where 


Assuming  a  linear  viscosity  temperature  relation 


=  c 


CO 


T 

7. 


and  a  unit  Prandtl  number,  the  system  of  differential  equations  governing  the 
boundary-layer  problems  in  plane  flow  is  reduced  to 


no 

=  e 
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=  0 


and  the  boundary  conditions  to 


=0  =  0 

at 

-1  =  0 
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=  CO 

#xere 
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H-N^ 

(2.1b) 
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In  arriving  at  these  equations ^  quantities  of  orders  and  Au,^/U  hi^es? 

than  the  remaining  terms  have  been  neglected#  Implicit  in  Eq«  (261)  is  'ttxe  fact 
that  Toithin  the  degree  of  approximation  afforded  by  the  boundary-layer  analysis^ 
the  normal  pressure  gradient  can  be  ignored#  For  the  neglect  of  the  pressure 
variation  across  the  boundary  layer ,  an  extremely  thin  boundary  layer  is  not 
strictly  necessary#  Rather,  the  neglection  depends  on  the  density  level  ■within 
the  hypersonic  boundary  layer  being  relatively  low#  Therefore  the  pressure  may 
be  taken  as  a  function  of  ^  alone  even  ■siien  the  boundary-layer  thickness  be- 

i(.«6 

comes  comparable  to  the  lateral  extent  of  the  outer  flow  # 

ETow,  consider  the  situation  vftien.  €  is  also  small#  Except  for  the  param¬ 
eter  e  itself,  all  terms  appearing  in  Eq*  (261)  will  remain  of  order  unity# 
That  is,  they  ■will  be  independent  of  e  in  the  limit  G  — >  0  #  In  this  limit , 
the  right-hand  member  of  Eq#  (26la),  ■ciiich  involves  the  pressure  gradient,  van¬ 
ishes  with  €  *  The  leading  approximation  to  the  solution  of  the  system 

Eq#  (2»1)  for  small  e  ■will  therefore  be  independent  of  e  and  is,  in  fact, 
governed  by  the  Blasius  equation 

Z  $  -f-  $  $  =0 

=  $'(o)  =  $ '  (co) -1^0 

with  ®  The  problem  of  determining  the  higher-order  approximation  is 

equivalent  to  fi.nding  a  correction  to  the  classical  flat-plate  solution  for  a 
sli^t  nonuniformity  of  the  outer  inviscid  flow#  The  nonhomo geneous  term  appear¬ 
ing  on  the  right  of  Eq#  (2#!)  reveals  clearly  the  role  of  the  wall-to-stagnation 
temperature  ratio  in  controlling  the  degree  of  approximation  of  local  boundary- 
layer  similarity#  Changing  Ty^ from  unity  to  zero  may  "typically  reduce  the 
correction  term  by  75^# 


Neglect  of  the  pressure-gradient  term  in  the  boundary-layer  momentum  equa- 

6 

tion  has  also  been  eiiq)loyed  by  Lees  to  siitgjliiy  the  problem  of  strong  shock 
boundary-layer  interaction,,  For  this  problem however  ^  a  self-siailar  solution 
exists  even  with  the  pressure  gradient  retainedo  The  theoretical  justification 
for  the  neglect  of  pressure  gradient  to  cases  involving  arbitrary  pressure  dis¬ 
tribution  on  the  basis  of  the  shock-layer  theories  has  only  been  pointed  out 

t  a* 

quite  recently  « 

The  mathanatical  problem  of  determiriing  the  higher-order  approximation  is 
one  that  involves  solution  to  a  linear  partial  differential  equation*  An  alter¬ 
native  -sitich  avoids  solving  the  partial  differential  equation  is,  of  course,  the 
momentum-integral  method  usually  adopted  in  boundary-lsyer  analyses.  For  the 
purpose  of  determining  pressure,  skin  friction,  and  surface  heat-transfer  rate 
for  the  present  problem,  this  method  may  be  sufficient,  as  has  been  demonstrated 

by  Stewartson''s  analysis  for  the  strong  interaction  problem  of  the  flat  plate  « 

h9 

The  scheme  arrived  at  recently  by  Moore  ,  ■which  accounts  for  the  first-order 
correction  to  the  local  similarity,  is  essentially  one  of  a  similar  nature. 

When  the  momentum-integral  method  is  en^loyed,  the  local  similarity  assumption 
is  no  longer  necessary,  even  as  a  leading  approximation  (for  small  6  ). 

In  the  leading  approximation  for  small  6  ,  the  surface  heat- transfer 

rate,  the  skin  friction,  and  the  displacement  thickness  of  the  boundary  l^er  are 
determined  explicitly  from  the  press\ire  distribution*  For  the  case  of  plane 
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where 


These  results j  in  view  of  the  approximations  involved 5  will  be  subject  to  errors 
of  order  €  and  A  U  o  The  usual  errors  inherent  in  the  boundary-layer 
equations  as  well  as  errors  from  possible  slipj  temperature- jump ^  and  vorticity 

effects,  will  presumably  be  no  larger  than  those  incurred  by  the  approximation 

13,l6,ltli  ^  , 

U.  .  The  magnitude  of  will  be  examined  later » 

Before  proceeding  to  the  analysis  of  the  inviscid  flew  region,  the  choice 
of  a  reference  temperature  for  the  determination  of  the  constant  in  the  linear 
viscosity-temperature  relation  will  be  discussed.  The  assumption  of  such  a  vis¬ 
cosity  law  may  seem  to  be  unrealistic  since  the  temperature  variation  within  the 
hypersonic  boundary  layer  is  quite  large  and  the  detail  of  viscosity  law  may  be 
important.  The  approximation  may  be  considerably  improved  by  choosing  the  con¬ 
stant  C  to  represent  correctly  the  value  of  viscosity  at  an  appropriate  refer¬ 
ence  temperature.  Instead  of  choosing  the  wall  temperature  as  a  reference,  as  was 
proposed  earlier  by  Chapman  and  Eubesin  and  adopted  later  by  others  for  hyper¬ 
sonic  applications,  we  recommend  a  reference  temperature  defined  by 


I  t  3 


(2oW 


whence 


This  choice  of  is  particularly  meaningful  when  local  boundary- layer  similar¬ 

ity  holds,  as  in  the  present  problem.  For  we  may  then  apply  Crocco*s  Integral 
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relation^  and  the  reference  temperature  7^  given  above  is  seen  to  be  simply  an 

+ 

average  temperature  across  the  hypersonic  boundary  layer,  ioe® 


1  “»  1 


The  skin  friction  and  the  heat-transfer  rate  on  a  flat  plate  under  uniform  pres¬ 
sure  based  on  of  Eq.  (2oU)  agree  exceedingly  well  with  results  based  on  the 
Young  formula'^  in  the  higher  Mach  number  range,  provided  the  Prandtl  number  is 
not  far  from  unity,  A  comparison  of  these  results,  including  the  results  based 
on  as  a  reference  temperature,  is  given  in  Table  I, 


An  examination  shows  that  7^  given  Eq,  (2,U)  turns  out  to  be  quite  close 
to  the  reference  temperature  proposed  by  Eckert^  , 
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SECTION  3 


LEADING-EDGE  BLUNTNESS  EFFECT  AND 
DETACHED  SHOCK  LAYER 

A  Model  for  the  Inviscid  Flow  Field 

The  necessity  of  aHo^Ming  the  shock  layer  to  detach  from  the  afterbody  of 

a  blunt  leading  edge  is  eridenced  by  the  breakdom  of  the  standard  shock-layer 

theory  around  the  shoulder  region  Tdiere  the  pressure  ■would  accordingly  fall  be- 

lo'w  zero  and  the  ^ock«layer  thickness  would  become  infinitee  The  essence  of 

otir  analysis  is  that  the  pressure  at  the  base  of  the  detached  layer  wiH  not  be 

39 

assumed  equal  -to  zero  5  a  prioid.^  as  ■was  done  for  a  free  layer  ,  "but  is  to  be 
determined  as  part  of  the  solution  to  the  problem* 

It  is  essential  to  realize  that^  wLthin  the  inviscid,  hypersonic-flo'w  re¬ 
gion  bounded  by  the  shock  and  the  body  (ieo*,  the  outer  edge  of  the  boundary 
layer)  downstream  of  a  small  blunt  nose,  there  is  an  inner  core  of  considerable 
extent  containing  gas  particles  -which  ha've  come  through  the  comparatively  forward, 
and  hence  stronger,  portion  of  the  shock  (see  Fig*  !)•  It  is  appropriate  to 
tem  this  region  an  entropy  layer,  since  its  specific  en-tropy  is  much  higher  than 
that  prevailing  near  the  adjacent  shock*  The  temperature  ■wLthin  this  layer  is 
accordingly  much  higher,  and  the  density  much  lower,  than  the  corresponding 
values  near  the  shock|  the  gas  within  the  entropy  layer  has  been  heated,  so  to 
speak,  by  the  strong  for-ward  portion  of  the  curved  shock*  It  is  natural  to  con¬ 
sider  the  remaining,  and  relatively  thin,  high-density  part  of  the  flow  field 
behind  the  shock  as  a  shock  l^er  in  the  usual  sense,  except  for  -the  fact  that  it 
is  no  longer  attached  to  a  rigid  body*  The  equivalent  body  surface  for  this 
shock  layer  over  the  relatively  slender,  downstream  portion  of  the  flow  field 
can  be  taken  to  coincide  with  the  outer  edge  of  the  entropy  layer,  (^ )  5 
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•which^  in  ■view  of  the  particle-isentropic  condition^  must  also  be  a  streamline 
(see  Fige  lb)«  Across  this  streamline  bo'undarya  the  pressure  must  be  contin~ 
uous» 


The  flow  structvire  -within  a  hypersonic  shock  l^er  has  been  treated  quite 
extensively  by  Freeman-^  and  others'^  "  «  Our  present  study  -will  deal  with  the 

entropy  layer,  the  outer  edge  of  which  is  the  base  of  the  shock  layer To  do  so 
we  need  a  knowledge  of  the  pressure  at  the  base  of  the  shock  layer  in  terms  of 
Y  »  The  leading  approximation  under  (b)  and  (c)  for  this  relation,  according 


to  shock" layer  analysis,  is  simply 


-p^  p  \J 


(3.1) 


Tiiich  represents  the  Ne-wtonian-plus-centrifugal  (Busemann)  pressure  formula  sim¬ 
plified  for  thin  bo  dies  <>  .  The  corresponding  formula  in  the  axisymmetric  case  is 


Restrictions  of  the  Present  Analysis 

In  view  of  the  high  temperature  attained  within  the  entropy  layer,  the 
local  velocity  U,  could  be  considerably  reduced  from  its  free-stream  value  U 
Since  the  assumption  of  ^  U  has  been  made  in  the  boun.dary-layer  analysis, 
and  will  be  used  later  in  the  treatment  of  the  entropy  layer,  an  assessment  of 
the  magnitude  of  is  important  and  will  be  made  prior  to  further  develop" 

mento 


It  may  be  noted  that  in  the  vicinity  of  the  bluht«nose  region,  the  inner  edge 
of  the  shock  layer  (vdrLch  in  this  instance  is  a  free  layer)  may  not  be  closely 
represented  by  a  streamline*  This  implies  that  streamlines  within  the  entropy 
layer  actually  come  from  various  points  along  the  curved  bow  shock  upstream 
(see  Fig*  la) ,  and  the  distribution  of  entropy  within  the  entropy  layer  is  not 
generally  uniform*  The  validity  of  the  assumption  that  the  inner  edge  of  the 
detached  shock  l^er  approaches  a  streamline  depends,  in  fact,  on  the  small- 
distance  assumption  (a) . 


For  the  purpose  of  examining  the  magnitude  of  in  the  entropy  layer  j 

we  consider  first  the  particular  case  of  a  flat  plate  under  the  influence  of  the 

strong  blast-wave  effect  produced  by  its  blunt  tip*  For  this  case,  the  surface 

pressure  distribution  is  quite  well  established  by  the  consistent  results  of  the 
19  20  2li 

blast-wave  theory  ^  ^  and  the  correlated  nximerical  calculations  based  on  the 

23 

characteristics  method  .  ¥e  can  therefore  determine  the  maximum  value  of  the 
perturbation  velocity  (for  arbitrary  /Kj  )  from  the  particle-isentropic  condition 
and  the  Bernoulli  relation.  Maximum  values  of  |  Aii/u\  for  /  *=  1,20,  1*1^0, 
and  1,667  are  presented  in  Fig,  2  as  a  function  of  n^/Jkt  *  "^he  results  indicate 
that  sufficiently  far  downstream  of  a  blunt  nose  (  ~  10) >  Au./u  remains 

sufficiently  small  to  permit  the  approximation  u.  U  . 

Since  in  most  cases. 


r*i 


where  Q  is  the  local  shock  angle  (which  gives  also  the  typical  local  flow 
angle) ,  a  conservative  estimate  for  the  perturbation  velocity  may  be  obtained 
from  the  Bernoulli  relation,  which  can  be  taken  as 


A  u. 


U 


A  y 


■vtiere  the  exponent  2  (y-  has  the  values  l/3,  h/l t  and  U/5  for  Y  -  1,20, 
1*U0,  and  1,667  respectively, 

¥e  may  now  state  more  explicitly,  in  teirnis  of  the  typical  flow  deflection 
angle,  the  requirements  >tiich  win  control  the  degree  of  accuracy  of  the  subse¬ 
quent  analysis  of  the  inviscid-flow  region,  as  well  as  of  the  boundary  layers 


At  Y/Jlit  -  10  >  \Au,/u\  ^  0,12,  O0I6,  and  0,26  for  7'  =  1,667,  Idi-O,  and  1,20, 
respectively*  In  the  calculation  for  Y  =1,20  ,  the  pressure  formula  given 
by  the  zero-order  theory  for  small  s  has  been  used. 
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^  Y 

(A)  must  be  negligible  in  comparison  with  unity 

(B)  Must  be  sufficiently  small  in  comparison  idth  unity 

(C)  6  must  -be  sufficiently  small  in  comparison  -with  unity 

Obviously^  condition  (a)  is  essential  for  the  assumption  iL=  U  at  the  outer  edge 
of  the  boundary  layer  and  mthin  the  entropy  layer  as  weHo  Condition  (B)  calls 
for  a  sufficiently  strong  shockj»  and  (c)  and  (B)  together  1)0.11  permit  treatment 
by  the  shock-layer  concepto  By  these  requirementss  the  admissible  ranges  of  the 
shock  (and  flow)  angle  9  ,  and  g  ,  are  necessaridy  restricteds  To  satisfy 
the  small-perturbation  assumption^  G  must  be  small^  but  not  so  small  that  the 
shock  cannot  remain  strongo  In  addition,  the  shock-layer  approximation  requires 

6  to  be  close  to  unity,  but  not  so  close  that  the  small-perturbation  assumption 
will  be  violated*  The  ranges  for  6  and  e  can,  of  course,  be  enlarged  by  in¬ 
creasing  M  and  by  including  (consistently)  terms  of  higher  orders  in  e  and 

Analysis  of  the  Entropy  Layer 

Our  analysis  will  be  based  on  the  fact  that  the  density  within  the  entropy 
layer  is  low  compared  to  that  in  the  free  stream,  "ice* 

essential  to  examine  the  con^tibillty  of  this  fact  with  the  requirements  set 
forth  above*  Since  the  density  depends  not  only  on  the  pressure  bub  also  on  the 
entropy,  w©  should  distinguish  two  situations*  In  one,  the  streamlines  over  most 
of,  or  at  least  a  considerable  part  of,  the  entropy  l^er  come  from  the  immediate 
vicinity  of  the  blunt  nose*  In  other  situations,  the  streamlines  from  the  blunt 
nose  region  reach  only  a  small  fraction  of  the  entropy  layer* 

The  rest  of  the  streamlines  in  the  entropy  layer  come  from  the  slender  but 
still  comparatively  forward  and  stronger  portion  of  the  curved  bow  shock* 
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In  the  first  situation,  the  specific  entropy  is  at  the  same  level  as  that 
prevailing  at  the  nose,  and  the  order  of  magnitude  of  the  density  may  therefore 
he  represented  for  particle-isentropic  flow  hy 


Comparing  this  with  requirement  (A),  we  see  that  the  local  density  is  indeed 

■55* 

small  in  comparison  with  the  free-stream  value  ■{dien  (A)  is  satisfied*  In  the 

second  situation,  the  density  distribution  given  ly  the  strong  blast-wave 

18,19,20,25,26 

theory  should  be  valid,  since  streamlines  over  most  parts  of  the 

entropy  layer  now  come  from  the  slender  portion  of  the  shockj  hence,  the  hyper¬ 
sonic  small-disturbance  theory  applies*  The  fact  that  the  blast-wave  solutions 
reveal  the  structure  of  a  detached  shock  layer  and  a  low- density  entropy  layer 
is  well  known*  The  ratio  this  situation  may,  in  fact,  be  represented 

by  (assuming  a  flat-plate  afterboclsr) 

P 


-  =  0 


e  ^ 


^ich  is  small  for  all  practical  purposes*  The  second  situation  may  perhaps  be 
realized  only  for  an  extremely  small  blunt  nose  (hence  large  1  )  together 
with  an  extremely  high  Mach  number*  In  either  situation,  the  supposition 

P 


4<  / 


(3.2) 


■within  the  entropy  layer  can  be  well  justified*  Its  error  is  presumably  no 
larger  than  that  inciirred  by  the  assiamption  u.  |9  ,  that  is ,  requirement  (A) . 

An  immediate  consequence  of  the  low  density  in  the  entropy  l^er  is  that 
the  variation  of  pressure  across  it  is  small  enough  to  be  negligible*  ¥e  can. 
examine  this  from  the  equation  of  motion  governing  the  momentum  component  in  the 


Note  that  under  requirement  (A) ,  6^^  «  €  for  e  4<  1 
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direction  transverse  to  the  layer.  The  fractional  variation  in  pressure  across 
the  entropy  layer  is  readily  seen  to  be 


=  0 

V  Y.  JJ 

=  0 


«  / 


(3.3) 


ce> 


and  is  therefore  negligible  according  to  our  requirements. 

In  view  of  the  above  conclusionj  we  can  now  employ  the  simplification 

u  ^  U 

r  (3oU) 

for  study  of  the  flew  pattern  within  the  entropy  layer.  The  equations  stating 
consei^mtions  of  mass  and  entropy  are 


U  ^  +  /V-  ^ 


9 


t 


(a  ^  -f  /r  ^ 


7 


Q/v' 


Y/ 


j 


can 


Eliminating  the  density  ^  and  making  use  of  the  simplification  on  u  and  ^  ^ 
these  two  equations  give  an  equation  for  ,  or  ,  This  equation 

be  integrated  to  yield 


r- 


(3.5) 


for  which  the  following  boundary  condition  at  the  inner  edge  of  the  entropy  layer 
(also  the  outer  edge  of  the  boundary  layer) 


At'/OL  ^ 


at  s  Yj^ 
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has  been  satisfied « 


Putting  d. Ict/^  for  >  “we  can  integrate  Eq«  (3«5) 


once  again  and  obtain  an  equation  for  the  streamline 


(3.6) 


where  rri  is  a  constant  differing  for  each  streamline »  This  equation  should 
hold  throughout  the  entropy  layer.  Applying  it  at  the  streamline  boundary 


-z  Yg  (')^)  s  'w®  arrive  at 

\yj^)  -  Yt  -  consUnt . 


(3.7) 


This  simple  integral  relation  between  the  pressure  and  the  cross-section  area  of 
the  entropy  layer  is  basic  to  our  subsequent  development.  The  corresponding 
integral  relation  in  the  axisyraraetric  case  is 
"  z 

Vfe  Yl  “P's  ~  constant  . 

b  J 

The  pressure  is  that  at  .the  base  of  the  shock  layer,  Y  ~  Ye  ® 
is  given  by  shock-layer  theory  once  is  specified.  This  relation  for 

from  shock- layer  theory  can  be  used  to  eliminate  in  Eq,  (3«7)^  and  we  can 

write  down,  at  least  in  principle,  an  ordinary  differential  equation  involving 
only  Yg  <■  The  shock- layer  theory  gives  the  Busemann  formula,  Eq,  (3.1),  as  a 
leading  approximation.  The  pressure-area  relation  consistent  with  this  approx¬ 
imation  is  therefore  reducible  to  (from  Eq,  (3.7)). 

(Y,  -  Yj(Y^Yj'/=  constant  . 

In  carrying  out  a  higher-order  analysis  (in  e  and  )  >  we  should  not 
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only  employ  a  pressure  formula  correct  to  the  order  desired  ,  but  also  work  con¬ 
sistently  with  a  more  accurate  pressure-area  relation  based  on  Eq.  (3<>7)» 

The  constant  on  the  right  of  Eq.  (3 <>7)  is  necessarily  related  to  the  drag 
of  the  blunt  nose,  as  a  consideration  of  the  balance  of  momentum  in  the  stream- 
wise  direction  will  readily  show.  With  the  help  of  the  equation  of  continuity 
and  making  use  of  the  small-perturbation  requirement  (A) ,  we  can  write  for  M  >  1 
the  equation  relating  drag  and  momentum  flux  as 


D, 


w 


-A  = 


■fO 


y-/ 


^14 


(3.8) 


where  is  the  drag  contributed  by  the  pressure  distribution  around  the  small 

blunt  nose’*l'  The  second  term  on  the  left  of  Eq.  (3o8)  then  represents  the  drag 
contributed  by  the  afterbody.  The  first  and  the  second  integrals  on  the  right  of 
Eq.  (3o8)  may  be  interpreted,  respectively,  as  the  internal  and  kinetic  energies 
associated  with  the  cross-flow  field.  As  a  result  of  the  small-perturbation  re¬ 
quirement,  the  streamwise  perturbation  velocity  u.  does  not  enter  into  the  above 
relation. 

In  view  of  the  drastic  difference  in  the  density  between  the  entropy  layer 
and  shock  layer,  it  is  convenient  to  write  the  second  integral  on  the  right  of 
Eq,  (3»8)  as 


The  first-order  e  correction  to  the  Busemann  formula  has  been  given  ex¬ 
plicitly  by  Freeman  in  Ref,  37, 

The  nose  drag  has  been  defined  as  the  total  drag  of  a  symmetrical  nose, 

n^^mely,  Y"  “  -  Yo  without  top-bottom  symmetiy,  P^/Z  should 

then  be  taken  as  the  contribution  of  one  part  of  the  surface  which  begins  at 
the  apex,  or  the  stagnation  point,  of  tte  nose. 
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‘^re 


The  lateral  extent  of  the  shock  layer  is  generally  of  an  order  higher  than  that 
of  the  entropy  layer,  i,e*, 

(K-KJ  =  -I 

whereas  the  ratio  of  density  in  the  entropy  layer  to  that  in  the  shock  layer  is 

As  a  result,  the  magnitude  of  the  first  integral  is  seen  to  be  smaller  than  that 
of  the  second  by  a  factor  j .  To  be  consistent  with  Eqs»  (3«.2)  and  (3»3)s 
we  may  omit  the  first  integral  in  comparison  with  the  second.,  Now,  if  only  zero 
and  first=order  terms  are  to  be  retained  in  the  momentum-drag  relation,  the 
kinetic-energy  integral  can,  of  course,  be  further  simplified,  and  Eqo  (3o8) 
may  be  written  as 


€  Of^  +  2  € 


r 


ry< 


\J‘ 


(3.8a) 


which  is  correct  to,  and  includes  terms  of,  the  first  order  of  (€  +  » 

The  afterbody  drag  integral  on  the  left  is  of  the  same  order  as  -p  ^  Since 

2  Z 

,  the  second  terms  on  both  sides  of  the  equation,  as  well  as 
—  ,  will  be  of  order  (’€  +  higher  than  the  first  integral 

on  the  right,  and  they  can  all  be  dropped  in  the  leading  approximation.  This 
agrees  with  the  observation  that  for  e  small  the  number  of  degrees  of  freedom 
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excited  becomes  very  large,  and  most  of  the  energy  gained  by  the  flow  field  ap¬ 
pears  as  internal  energy,,  which  is  the  remaining  term  on  the  right. 


Since  the  pressures  in  the  shock  and  the  entropy  layers  belong  to  the  same 
order,  we  can  replace  the  location  of  the  shock  by  the  base  of  the  shock  layer 
in  the  zero-order  approximation,  and  the  momentum  relation  Eqo  (3o8),  to  this 
approximation,  is  reduced  therefore  to 

That  is, 

(X-Yj(Y»r;)'=  € 


This  is  the  pressure-area  relation  pertaining  to  the  zero-order  theory  previously 
given,  but  with  the  constant  now  determined.  For  the  axisymmetric  case,  the 
corresponding  relation  is 

n 

Z 


y/-  yJ]  [r;  V  Y.  Yj/z]  =  2  .  u 


For  the  zero-order  approximation,  it  seems  that  consideration  of  the 
momentum- integral  relation  not  only  determines  the  constant  of  Eq.  (30?)^  but 
also  furnishes  an  alternative,  and  apparently  simpler,  derivation  of  the  pressure- 
area  relation.  For  the  higher  approximation,  however,  we  may  find  it  convenient 
to  use  Eq.  (3o8)  strictly  as  a  subsidiary  condition.  The  constant  in  the 
pressure-area  relation,  Eq.  (3<.7)  j  can  be  determined  by  passing  to  the  limit 
in  Eq.  (3o8). 

¥e  may  note  that,  on  the  basis  of  Eq.  (3e9)>  the  pressure  ,  i.e., 

(Ye  YJ)  ,  can  never  become  negative.  ThisY is  to  be  expected  since  the  shock 
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layer  in  our  treatment  is  no  longer  required  to  follow  the  body  surface.  However ^ 
Eq.  (3.9) s  as  well  as  Eq,  (3.7),  has  not  excluded  the  possibility  of  an  attached 
shock' layer,  and  also  contains  the  free  layer  as  a  limiting  solution.  ¥e  can  see 
this  readily  by  examining  Eq.  (3.9)  for  the  nearly  limiting  situation  of  an  ex¬ 
tremely  small  nose.  The  constant  on  the  right  of  Eq,  (3.9)  in  this  case  then  re¬ 
duces  to  a  very  small  but  positive  number.  Now,  for  -p>^  to  be  finite,  it  is 
clear  that  the  shock  layer  must  remain  essentially  attached  to  the  body,  i.e,, 

Y  <5^  Y  •  soon  as  approaches  zero,  the  entropy  layer  will,  accordingly, 

begin  to  thicken  and  the  shock  layer  will  detach  to  form  a  free  layer. 

Discussion 

On  the  basis  of  the  foregoing  analysis,  the  fields  of  pressure  and  velocity, 
as  well  as  the  locations  of  the  shock  can  all  be  consistently  determined  under 
requirements  (A) -(C),  ¥e  have  arrived  at  these  results  without  the  need  of 
simultaneously  knowing  the  density  distribution  within  the  entropy  layer,  except 
for  the  fact  that  the  density  level  of  the  entropy  layer  is  comparatively  low. 

On  the  other  hand,  the  present  theory  actually  affords  no  solution  for  the 
entropy  fieldj  to  find  the  density  or  entropy  in  the  entropy  layer,  certain  de¬ 
tailed  knowledge  of  the  flow  field  around  the  blunt-nose  region  would  be  neces¬ 
sary,  It  imist,  however,  be  pointed  out  that  the  determination  of  density  may  no 
longer  be  a  trivial  matter  for  the  present  problem  if  the  assumption  of  constant 
specific  heats  were  removed.  For,  under  this  more  general  situation,  the  value 
of  density  would  be  essential  for  determining  the  local  thermodynamic  properties, 
hence  the  local  value  of  £■  ,  which  in  turn  will  determine  the  pressure  and 

velocity  fields.  In  this  case,  it  will  be  more  appropriate  to  redefine  %  as 
6  =  •^/2/?^and  the  relation  corresponding  to  Eq,  (3.9)  becomes  (V  ®  0,  or  1)  » 
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The  results  of  the  above  analysis  confirm;,  in  a  certain  senssj  the  validity 

l8  20 

of  the  blast-wave  theory  for  the  tip-bluntness  problem"^  “  o  In  spite  of  the 
fact  that  the  blast-wave  theory  of  Refs,  18,  19  and  20  does  not  take  proper  ac¬ 
count  of  the  presence  of  the  entropy  layer,  and  gives  erroneous  values  for  the 
entropy  and  density  near  the  body  surface,  the  pressure  and  velocity  determined 
by  the  present  analysis  of  the  entropy  layer  may  be  identified  as  those  of  the 
blast-wave  theory,  as  an  inspection  of  the  governing  Eqs,  (3 <>7)  and  (3»8)  will 
readily  reveal. 

Our  general  approach  to  the  inviscid  leading-edge  problem  is  also  similar 

21  22 

to  a  treatment  given  previously  by  CheiVQri  ^  ,  as  pointed  out  earlier.  In 

Cherr^ri’s  treatment,  which  is  essentially  a  momentum- integral  method,  use  is 
made  of  U  and  The  basic  equations  employed,  written  here  for 

steady  plane  flows,  are 

D, 


(3.10) 


We  note  that  the  first  of  Eq,  (3»10),  except  for  the  first  term  on  the  right-hand 
side,  is  almost  the  same  as  Eq,  (3«8a).  It  is  evident  that,  as  far  as  the  lead¬ 
ing  approximation  for  small  g  and  / y/ 0^  is  concerned,  the  above  equation 
and  Eq,  (3.9)  obtained  by  the  present  study  are  essentially  equivalent,  When 
the  higher-order  approximations  are  considered,  we  must  compare  Chernyi's  equa¬ 
tion,  i,e,,  Eq,  (3,10),  with  our  Eq,  (3,7),  We  may  note  that,  implicitly  in 
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Eq.  (3«10)j  the  outer  edge  of  the  entropy  layer  has  been  replaced  by  the  shock 
itself,  and  the  Busemann  pressure  formula  has  been  used.  Hence,  it  is  doubtful 
that  Eqo  (3ol0)  would  yield  a  consistent  correction  to  the  leading  approximation.. 
At  any  rate,  the  two  treatments  are  seen  to  be  quite  different  when  higher-order 
approximation  is  considered. 

The  entropy  layer  downstream  of  a  blunt  nose  has  also  been  studied  re¬ 
cently  by  Guiraud^^^^^,  By  an  ingenious  mathematical  technique  on  which  the 
writer  is  not  in  the  position  to  elaborate,  Guiraud  has  arrived  at  an  equation 
which  can  be  readily  identified  as  the  pressure-area  relation,  Eq,  (3«7) ^  given 
above.  However,  Guiraud  does  not  use  the  momentum-drag  relation,  Eq.  (3<.8)  or 
Eq,  (3.8a),  to  determine  the  constant  term  on  the  right  of  Eq.  (3.7).  Instead, 
he  relates  the  term  to  an  integral  involving  shock  curvature  in  the  nose  region. 
Hence,  this  term  cannot  be  explicitly  determined  until  the  complete  problem  in¬ 
volving  the  nose  region  is  solved.  Strictly  speaking,  therefore,  the  results  of 
Guiraud 's  study  in  Hefs.  ii2  and  1x3  do  not  fall  into  the  same  class  as  the  blast- 
wave  theory. 

The  reason  for  the  departure  of  Guiraud’ s  treatment  from  the  present  one  at 
the  final  stage  of  the  analysis  may  be  related  to  the  remarks  given  in  Ref,  ii2s 
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tftiereas  the  blast-wave  theory  provides  .  ^  oC  «  It  was  thus  concluded  that 
the  momentum-drag  integral  relation  vhich  leads  to  the  blast-wave  result  cannot 
be  used  to  determine  the  constant  since  it  appears  to  contradict  the  pressure- 
area  relation^  Eq<>  (3o7)^  xmless  / — ^1  «  It  must  be  notedj  however^  that  no 
streamline  in  a  self-similar  hypersonic  flow  field,  except  for  those  on  the  body, 
can  be  similar  to  the  shock  surface,  i.e®,  Y1  cC  Y  cannot  be  true.  Otherwise, 
there  could  be  no  flow  into  the  region  between  the  shock  and  the  body.  In  fact, 
Eq.  (3»7)  states  just  this  fact.  For,  according  to  the  blast-wave  results, 

Yg  oC  and  oc  ^  3  "the  streamline  boundary  for  the 

entropy  layer  will  then  be  given  by  the  pressure-area  relation  as 

.  ,  -///  Vi/  /  Z/s] 

Y  oc  oc:  J 


The  consistency  between  Y,  Y  small  but  nonvanishing 

e  can  indeed  be  checked  by  the  results  of  the  higher-order  analysis  of  the  shock- 
layer  theory^"^. 

Implicit  in  the  foregoing  discussion  is  the  assumption  that  the  physically 

meaningful  solution  to  Eq.  (3o7),  vhen  the  leading-edge  effect  predominates,  is 

the  blast- wave  solution  V  oc  «  The  fact  that  the  blast-wave  solution 

's 

does  not  necessarily  represent  the  complete  solution  to  Eq.  (3o7)  can  be  seen  by 
examining  the  corresponding  equation  for  small  e  in  the  case  of  the  flat-plate 
afterbody,  that  is, 

f 

2 
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The  complete  solution  of  this  equation  can  be  readily  found»*  However^  the  only 

particular  solution  ivhich  is  compatible  isith  the  small" disturbance  requirement 

turns  out  to  be  V  oC.  *  A  detailed  discussion  of  this  point  is  given  in 

Appendix  1*  For  reasons  of  consistency,  the  solution  to  Eqso  (3*7)  or  (3*9) 

iNhen  bluntness  dominates  Td-ll  be  taken  to  be  throughout  the  present 

's 

analysis • 

In  Tiihat  follows,  we  shall  apply  the  results  of  the  above  analysis  to  the 
problem  involving  both  tip  bluntness  and  boundary-layer  displacement  effects® 

%  9 

The  zero-order  theory  for  small  €  and  large  M  Q  will  be  developed  in  de¬ 
tail® 


*  The  existence  of 
attention  by  Ifr® 


a  complete  solution,  in  this  case,  was  called  to  the  writer's 
J«  P«  Guiraud  in  a  recent  communication® 
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SECTION 


THE  ZERO-ORDER  THEORY  FOR  INTERACTION  OF  THE  BOUNDARY  LAYER 
AND  ITS  OUTER  FLOW  UNDER  THE  INFLUENCE  OF  LEADING-EDGE  BLUNTNESS 


The  Governing  Equations  for  Small  e  and  Large 

For  small  €  }  the  boundary  layer  has  been  shown  to  be  governed  by  local 

similarity and  assuming j  in  addition,  the  shock  wave  involved  is  strong,  the 
simple  pressure-area  relation  Eq,  (3o9)  from  the  analysis  of  the  entropy  layer 
should  also  hold.  Combining  these  results,  the  leading-edge  problem  of  hyper¬ 
sonic  flows  subject  to  the  tip  bluntness  and  boundary- layer  displacement  effects 
can  be  treated  in  a  rather  straightforward  fashion. 

To  apply  Eq,  (3o9)  to  this  problem,  the  inner  edge  of  the  entropy  layer, 
f  will  be  taken  as  an  equivalent  body  surface  which  should  account  for 
the  displacement  effect  of  the  hypersonic  boundary  layer.  Inasmuch  as  the  assump- 
tion  LCj  U  is  valid  ^  ,  the  outer  edge  of  the  boundary  layer  may  be  approx¬ 
imately  represented  by  the  displacement  thickness  (5*^  because  of  the  extreme 
low  density  resulting  from  viscous  heating.  In  this  instance,  we  can  replace  Y 

^  ^  local  similarity,  <5'^  is  given  by  Eq.  (2,3)  explicitly 

in  terms  of  the  pressure  distribution,  which  is  now  provided  by  the  Busemann 
formula . 


The  equation  governing  the  zero-order  approximation,  that  is,  the  leading 


approximation  for  small  e  and  large  M ^  j  is,  for  the  plane  flow. 


[Ye-  Y.  -  (Ye  t'J  -  ^2 

6^  ^  €  (0,  664  +  qs 

^  ^  lo/  (Yer.’j' 


(U.i) 
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After  solving  Eq.  (l|ol)  for  s  we  can  determine  the  pressure  from  the 

Busemann  formula  and  the  heat=transfer  rate  as  well  as  the  skin  friction  from 
Eqo  (2  <,3)0  The  rudimentary  manner  in  which  the  three  flow  regions  ■=<=  boundary- 
layer  ^  entropy  layer.,  and  shock  layer  ==  interact  is  contained  in  the  ordinary- 
differential  equation  given  above.  In  the  absence  of  tip  bluntness,  we  have  only 
the  shock  layer  and  the  boundary  layer.  The  first  of  the  above  equations  then 
becomes 

and  Eq,  (U,l)  reduces  to 

Near  the  origin^  /)L  -=.  +  0  ,  Eq.  (U.l)  admits  the  singular  solution  » 

This  is  the  correct  behavior  when  the  blast-wave  effect  of  the  leading  edge  is 
19  20 

strong  ^  j,  and  confirms  that  near  the  leading  edge  the  displacement  effect  of 
the  boundary  layer  gives  way  to  the  more  powerful  effect  of  leading-edge  blunt¬ 
ness,  For  cases  which  do  not  involve  bluntness,  Eq.  (li.la)  reveals  singular 
behavior  of  another  types  as  — >0  $  Y  ^  ^  »  This  is  the  behavior  to  be 

expected  when  the  displacement  effect  predominates-^^  ^  .  With  these  singular 
solutions  as  leading  terms,  solutions  to  Eq,  (l|,l)  can  be  obtained  by  a  forward 
integration  from  the  origin.  These  asynptotlc  solutions  for  Y-  5  as  well  as 
the  pressure  and  the  heat  transfer,  will  be  examined  subsequently  for  flat-plate 
afterbodies. 

For  Eq,  (ii-.l),  we  shall  assume  that  as  <7^ — >0  }  ,  where  cr  >  2/3  , 

For  Eq,  (l|,la) ,  we  ass-ume  that  as  Ti' - ,  Y^  s  where  cr  >  ,  In 

a^  case,  assuming  a  regular  afterbody  shape,  i,e,,  a/  <9^  as  , 

will  be  sufficient. 
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As  is  obvious  from  the  requirements  {k)~{C),  the  error  in  the  zero-order 


theory  is  presumably  of  the  order 

0 


Y-l 


^  T 


If  terms  of  the  next  orders  are  consistently  included  in  the  analyses  of  the  en¬ 
tropy  layer  and  of  the  boundary  layer^  in  the  manner  previously  described,  the 
first  two  error  terms  may  presumably  reduce  to  and  • 

In  view  of  their  frequent  occurrence  in  the  later  development,  it  is  con¬ 


venient  to  introduce  two  new  variables  for  .the  plane  flow, 

i 


K.  s 


(lt.2) 


where 


Re  — 


Except  for  the  factors  6  and  e  (0«66J4  +  1.73  %i/Tq  )  j  these  two  variables  are 

l-2ii 

those  employed  previously  for  the  bluntness  and  the  displacement  problems  , 

In  what  follows,  the  zero-order  theory  will  be  applied  to  the  problem  of  a 
flat  plate  in  hypersonic  flow  subject  to  the  three  effects  of  bluntness,  dis¬ 
placement,  and  angle  of  attack.  In  order  to  see  clearly  the  manner  in  which  these 
effects  will  act  in  combination,  the  simpler  problems  involving  only  two  simul¬ 
taneous  effects  (nanielys  ^  ^  “■  Mcxl  j  K  ~  Mo<.)  will  first  be  singled 

out  for  discussion. 

The  Flat  HLate  with  Bluntness  and  Displacement 

For  this  piublem,  Eq,  (^,1)  can  be  reduced  to 


30 


Z  C2?'/-/zh'  =/ 


(U.3) 


where  the  new  variables  are 

Z  -  8M 

c  =  /<^ 

Note  that  both  (- sind 


(0„66i;  +  1„73  M^C/ARe 


are  independent  of  o  Therefore  j  H  and  ^  are  proportionals,  respectively ^ 


to  and  /p  0  Let  (j)  and  A  =  cL4>/d^  o  Equation  (ii.<,3)  may  then  be 

A  -  fX  =  / 


written 


which  can  be  Integrated  by  a  separation  of  variables  to  give 

0  =4-  [/T  -  x/z  -  A,'  T/t-  t^T)] 


Vs 


This  result  provides  the  desired  behavior  at  the  origin^  namely 

o  To  obtain  a  relation  between  ^  and  (j)  y  "we  note  that 

A 


/  ^  ^  M  J  / 

CL<^  =  •  ,  J-  d,  <p  - 


d  (/> 


A 


and  the  solution  of  Eq<,  {h>3)  can  finally  be  brought  to  the  parametric  form 

z  =  2  /T  ^  /X  —  ^z  X  YX ^3  ~  yZfv  (7  +  yx”  ^ 


^  ^  ^  /r;%  9  ^  /r) 


9 


All  A) 
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This  result  for  s  which  may  also  be  taken  as  the  shock  location 


%/3 


The  corre- 


ie  presented  in  Fig,  3(a)  as  ^e/^e 

spending  surface  pressure  and  the'  surface  heat-transfer  rate  are  given  in  (b) 
and  (c)  of  the  same  figure  as  ^  ®  (o,66ii  +  1,73 

¥e  can  see  that  y^  varies  as  as  — >■<»  , 

and  as  as  -^-—>0  ,  These  two  limits  correspond  to  situations  in  which 

displacement  and  bluntness  dominate^  respectively,  agreeing  with  strong-interac¬ 
tion  theory  in  the  first  case  and  strong  blast-wave  theory  in  the  second.  One 
important  feature  here  is  that-  all  parameters  and  spatial  variables  are  grouped 
together  to  form  one  single  variable  0  From  Fig,  3  we  see  that 

transition  from  the  strong  blast-wave  to  the  strong-interaction  regimes  occurs 
in  the  range 

Id''  <  zVk/'  <  2 


A  variable  analogous  to  ,  which  controls  the  combined  effect  of 

blunt ness  and  displacement,  has  been  previously  pointed  out  by  Lees  and  Kubota 

2 

in  a  discussion  assuming  a  weak  nose  effect  , 

The  appearance  of  and  in  the  nondimensional  forms  of  pressure 

and  heat-transfer  rate  introduced  above  result  from  the  use  of  the  Busemann  for¬ 
mula 


Retaining  the  factor  Y  ,  instead  of  replacing  it  by  unity,  is  in  keeping 

'  <o 

with  the  convention  of  existing  shock- layer  theories^^”^^.  Furthermore,  it  per¬ 
mits  €  to  represent  conditions  behind  the  shock  even  when  is  not  small. 
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It  is  of  interest  to  examine  more  closely  the  asymptotic  solutions.  They 


,  respectively ,  as 


can  be  written,  for  large  and  small  >  respective! 

MYj/*;  =I.ZI)%  >  l-ilKY’ 

X^^O.ZIS/^  >  O.II9(K.f 


(lt.5) 


The  appearance  of  Y  in .the  above  strong- interaction  solutions  indicates 
clearly  the  roles  of  €  and  the  wall  temperature  level.  In  view  of  the  factor 
(0,661;  +  1.73  T^/r^  )  ±n  Xg  s  the  boundary- layer  displacement  effect  will  de¬ 
pend  strongly  on  the  wall  temperature.  An  interesting  observation  is  that  if 
we  simply  add  the  two  asymptotic  solutions  for  s  the  sum  differs  from 

the  exact  solution  in  the  transition  region  by  less  than  12^.  This  seems  to 

lend  support  to  the  empirical  linear- combination  laws  proposed  previously  by  cer- 
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tain  investigators  ,  at  least  when  €  is  small. 

¥e  may  now  compare  the  zero-order  solution  in  the  two  limits  with  the  cor¬ 
responding  self-similar  solutions  which  do  not  assume  close  to  unity.  In 
Tables  II  and  III,  results  from  a  few  of  these  existing  solutions^^^'^*^^^^^^^^ 
for  'Z'  =  l.UO  are  tabulated  and  compared  with  those  of  Eq.  (ho^) »  Some  of  these 
results  are  also  shown  in  Fig.  3.  conclude  from  these  comparisons  that 

as  an  approximation  for  -Wie  surface__pres^re  aM  skin  fraction  when  Z  ~  ? 

the  present  zero-order  solution  may  not  be  adequate ,  and  a  higher-order  solution 
(in  €  )  appears  desirable.  For  the  surface  heat-transfer  rate  as  well  as  the 

shock  shape,  however,  the  zero-order  theory  works  reasonably  well,  in  spite  of 
the  fact  that  Z  in  this  case  is  actually  not  too  close  to  unity. 
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The  Inclined  ^Flat  Plate  with  Boundary°Layer  Displacement 

With  bluntness  excluded,  Eq.  (liola)  can  be  transformed  to 


(e  7  =  ! 


{k.k) 


where  the  variable  2  and  ^  stand  for 

2 


J.  4  ^ 

z  ^  4M  OC 


5  = 


The  angle  of  attack  of  the  plate,  oc  ,  is  positive  for  compressive  flow  over 

the  plate o  The  negative  sign  in  front  of  C  EQo  (li»6)  applies  to  oc  >  0  » 

In  view  of  the  strong  shock  requirement ,  our  solution  for  o<^  <1  0  pertains  only 

to  the  situation  when  the  boundary  layer  is  thick  enough  to  induce  a  strong  bow 

shock.  This  implies  that  the  region  of  validity  is  confined  to  certain  upstream 

portions  of  the  flow  field,  and  to  a  rather  small  range  of  negative  oC  ,  .  Again, 

there  is  only  one  variable  governing  the  problem.  In  this  case,  it  is 

Equation  (1^,6)  appears  to  be  quite  different  from  that  underlying  the  analysis 
1  7  10 

of  Pai  and  Shen  ^  ^  ,  The  difference  presvimably  arises  from  the  use  of  a  dif¬ 

ferent  formula  for  the  pressure. 

For  the  case  of  c<  >  0  the  asymptotic  solution  of  Eq,  (1|,6)  for  large 

,  i,e,,  large  ^  ,  is  simply  that  for  a  wedge  without  displacement 

effect,  ioe, 

I 

(h.7) 


e(o. (,64-1.73^)^,//^  7^  =  oc/m^ 


3h 


For  small  argument  of  the  variable,  the  strong -interaction  solution  results,  as 
given  previously  in  Eq.  When  the  wedge  angle  is  sufficiently  small,  or  ^ 

small,  the  incidence  effect  amounts  to  only  a  small  perturbation  on  this  solu¬ 
tion,  ¥e  may  then  write  Y  or  y  as 

£  s 

(1-8) 


This  development  near  the  sharp  leading  edge  is  valid  for  both  od  >  0  and 
oC<0  ,  and  reveals  a  pattern  of  asymmetry.  The  curves  for  Cd  >  0  presented 

,  -it 

in  Fig,  4  give  a  smooth  transition  from  the  regime  of  strong  interaction  to  that 
of  a  usual  hypersonic  wedge  flow,  Cuives  for  od  <.  0  in  the  same  figure  show  a 
gradual  deviation  from  the  Stewartson-Lees  solution  as  we  go  downstream.  For  the 
shock  shape  and  the  pressure,  the  range  of  transition  for  od  >  0  appears  to 
occur  in 

M^ccy^z^  <  10 


whereas  for  the  heat-transfer  rate,  the  value,  corresponding  to  the  pure  wedge- 

^  z/ 

flow  solution  is  closely  approached  as  early  as  M  “  2  ,  From  the  agree¬ 

ment  of  our  zero-order  solutions  with  the  self-similar  solution  for  Y  =  loliO  in 
the  two  limits  (see  Fig,  U),  we  may  anticipate  the  zero-order  results  to  be  rea¬ 
sonable  for  heat  transfer,  pressure,  and  shock  shape,  provided  the  surface  tem¬ 
perature  is  substantially  lower  than  the  stagnation  temperature. 


The  numerical  results  of  Fig,  i;  were  obtained  by  digital-machine  computation 
(IBM  >  using  an  interval  of  100  points  per  unit  of 


35 


The  Inclined  Flat  Plate  with  Tip  Bluntness  and  the  Corresponding  Problem 
of  a  Blunted  Cone 

For  the  inviscid  problem  of  the  blunted  flat  plate  at  incidence  with  neg¬ 
ligible  boundary-layer  displacement  effect,  we  have  from  Eq.  (I4.0I)  the  simple 
differential  equation 

(Z  +  C)(Z  ZO'  =  /  (1,.9) 

where  the  variables  2  and  ^  now  stand  for 

At 
At 

In  this  case,  we  again  use  the  same  convention  for  oc  . 

For  the  purpose  of  later  discussion,  we  will  examine  in  detail  the  asymp¬ 
totic  behavior  of  the  solution  of  this  inviscid  problem  for  small  and  large  ^  » 
Near  the  leading  edge,  we  find  from  Eqo  (i|-o9)  the  desirable  result 

which  holds  for  both  positive  and  negative  and  reveals  an  asymmetrical  per¬ 
turbation  of  the  blast-wave  solution.  At  the  far  downstream  end,  where  the 
effect  of  the  surface  incidence  dominates,  we  will  develop  the  asymptotic  solu¬ 
tion  only  for  oC  y-  0  „  Equation  (U»9)  then  yields  for  large  ^ 

z  ^  a.u) 

where  /\  and  Q  are  constants,  and  and  are  Bessel  functions  of  the 

first  and  second  kind  of  certain  order  (T  «  In  view  of  the  first  and  second 
terms  we  see  that  the  shock  will  tend  to  the  wedge  surface,  but  never  quite  reach 
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it,  as  we  go  downstream.  Most  interesting  is  the  oscillatory  behavior  of  the 
solution.  According  to  the  asymptotic  behavior  of  the  Bessel  functions,  the 
third  group  of  terms  of  Eq.  (li.lO)  is  equivalent  to 

*B) 

The  corresponding  behavior  of  the  pressure  is  then 

Evident  from  Eqs.  (It-.ll)  and  (l4.,12),  this  oscillatory  decay  of  the  tip-bluntness 
effect  will  be  more  readily  revealed  in  the  pressure  distribution  than  in  the 
shock  shape.  The  results  obtained  from  numerical  integration  of  Eq,  {ho9)  by 
digital  machine  computation'*"  are  presented  (as  the  full  curves)  in  Fig,  The 
oscillatory  decay  of  the  tlp=bluntness  effect  is  indeed  exhibited  by  the  pressure 
distribution  for  OC  >  0  , 

It  may  be  of  interest  to  point  out,  in  this  connection,  that  the  solution 
to  Eq,  can  be  approximated  surprisingly  well  by  the  Equation 

=(az'/=/  +  7^[A59  f2/y;-a7/ r^(2/r^ 

which  is  also  included  (as  dotted  curves)  in  Fig,  5-  The 'function  "2  given 
above  not  only  agrees  with  the  asymptotic  solution  Eq.  (ii.ll)  for  largb  ^  , 

but  also  has  a  singularity  at  ^=0  identifiable  as  that  of  the  blast-wave  solu¬ 
tion,  In  addition,  it  satisfies  an  integral  relation  of  the  differential  equa- 

An  interval  of  100  points  per  unit  of  (C)  was  used  in  the  computation. 

The  heat-transfer  rate,  which  is  not  presented,  may  be  obtained  in  terms  of 
the  pressure  from  Eq,  (2,3). 
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tion  Eq<.  naxnely' 


r 

r '  -/] 

i 

=/ 


The  physical  significance  of  the  oscillatory  decay  observed  from  the  above 
analysis  should j,  of  course ^  be  taken  with  cautiono  The  result  of  the  zero-=order 
solution  in  this  respect  is  evidently  connected  with  the  passage  to  the  mathe¬ 
matical  limit 5  >1  «  Apparently  the  oscillation  represents  a  repeated  re¬ 

flection  or  skipping  of  the  reattaching  shock  layer  as  it  reaches  the  vicinity 
of  the  surface  of  a  very  slender  afterbody  at  infinite  Mach  number.,  The  phenom¬ 
enon  may  hence  be  interpreted  as  a  form  of  glancing  reimpingement  of  the  shock 
layer.  From  the  very  fact  that  this  oscillation  is  a  decay  phenomenon,  its  de¬ 
tail  may  probably  be  substantially  altered  by  higher-order  6  corrections.  For 
loUOj  in  reality  we  may  not  have  ary  oscillation  in  the  pressure  decay,  since 
the  amplitude  of  the  fluctuation  shown  is  never  large. 

For  the  purpose  of  comparing  the  afterbody  effects  in  two  and  three  dimen¬ 
sions,  the  results  of  a  similar  analysis  applied  to  a  blunted  slender  cone  will 
be  briefly  discussed.  The  basic  differential  equation  corresponding  to  Eq,  (U.?) 
for  this  axisymmetric  problem  is 

(z”- z  z'/z )  =  / 


with 


z 


t 

t 


This  approximate  analytic  representation  of  the  solution  to  Eq,  (U,?)  was 
pointed  out  to  the  author  by  FTofessor  H.  Pollard  of  Cornell  University, 
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where  oc  now  stands  for  the  half-cone  angle,  t  the  nose  diameter,  and  Jk/  the 
nose  drag  coefficient 


i  A  U"irt/i 


The  peculiar  manner  in  which  the  shock  layer  approaches  the  afterbody  is  more 
markedly  revealed  in  this  case  than  in  the  case  of  the  blunted  wedge.  From  the 
resulte  for  shock  shape  presented  in  Fig,  6(a) ,  we  see  that  the  changeover  in 

^  = /  occurs  rapidly,  as  if  the  detached  shock  layer 
remains  essentially  unaware  of  the  presence  of  the  afterbody  until  reaching  the 
immediate  neighborhood  of  the  cone  surface.  On  the  other  hand,  the  corresponding 
surface  pressure,  as  well  as  the  surface  heat-transfer  rate  (calculated  on  the 
basis  of  local  similarity  for  the  laminar  boundary  layer) ,  given  in  (b)  and  (c) 
of  the  same  figure  show  a  much  earlier  influence  of  the  afterbody.  The  oscilla¬ 
tion  in  pressure,  decaying  as  a  Bessel  function,  is  again  found,  but  in  a  rather 
manifested  way.  We  should  note  that  in  this  axisyrametric  case,  the  first  pres¬ 
sure  undershoot,  and  the  oscillatory  cycle  immediately  following  it  are  so  pro¬ 
nounced  that  their  occurrence  may  not  be  too  surprising  even  when  £  is  not  as 
small,  and  the  shock  not  as  strong,  as  required  by  the  zero-order  approximation. 
This  peculiar  behavior  of  the  pressure  distribution  may  be  quite  meaningful  for 
studies  of  the  influence  of  bluntness  on  boundary-layer  instability.  The  pres¬ 
sure  undershoot  has,  in  fact,  been  observed  experimentally  for  slender  blunted 
cones  at  /vj  ”  6.85  in  ain 

^  22 

The  undershoot  of  cone  pressure  has  been  predicted  previously  by  Cherryi  ^  , 

using  an  integral  method  discussed  in  Section  2.  For  the  inviscid  flow  past  a 
blunted  cone  with  '/  =  l.ltO,  Cherryi  finds  a  strong  overshoot  following  the  pres- 


the  vicinity  of  oC 


i/a 
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sure  undershoot.  The  Bessel  function  decay  is,  nevertheless,  not  found.  In  a 

21 

similar  treatment  for  the  blunted  wedge  ,  there  is  no  pressure  undershoot,  ac¬ 
cording  to  Chernyi’s  result.  The  difference  betvreen  the  present  solutions  and 
Chernyi's  presumably  results  from  the  difference  between  Eqs.  (3«9)  and  (3<>10). 

The  Combined  Effect  of  Bluntness,  Displacement  and  Body  Thickness 

For  the  flat-plate  problem  involving  all  three  effects  (Z  > 

€.  ® 

,  we  find  that  a  new  parameter  appears  in  the  reduced  form  of  the  solution. 
If  we  choose  to  use,  for  example, 

Eq.  (U.l)  for  this  problem  becomes 

(z.  -  -/  (U.13) 

where  p  is  a  parameter  (independent  of  )  given,  for  the  wedge  flows,  as 

n  — 

Figure  7  gives  the  results  of  Integration  for  the  surface  heat-transfer  rate  in 
the  same  variables  as  used  for  the  zero-incidence  case,  i.e..  Fig.  3(c).  Shown 
are  curves  corresponding  to  five  positive  values  of  p  ,  namely,  P  =  0,  10  , 

10  1,  and  10.  The  curve  for  P  ==  0  has  been  previously  given  in  Fig.  3(c). 

Similar  calculations  may  also  be  carried  out  for  negative  values  of  P  , 

From  these  results  we  can  see  the  role  of  the  wedge  angle  in  modifying  and 
controlling  the  combined  effect  of  the  bluntness  and  boundary- layer  displacement. 
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The  slight  oscillatory  mode  revealed  in  the  curve  at  the  upper  left  corner  of 
Figo  7  for  P  =  10  corresponds  to  the  oscillatory  pressure  behavior  discussed 
previously  for  the  inviscid  wedge  problem.  The  results  show  that  when  P  is 
not  too  larges  the  oscillatory  phenomenon  will  subside  as  a  result  of  the 
inviscid-viscous  interaction. 

Hypersonic  Similitude  Eased  on  the  Zero-Order  Theory 

From  the  flat-plate  problems  we  see  that  the  zero-order  solution  is  redu¬ 
cible  to  a  function  involving  only  om  variable  and  the  parameter  P  ,  Viscous 
hypersonic  flows  over  a  blunted  inclined  plate  will  accordingly  be  similar  for 
the  same  P  „  For  thin  bodies  other  than  flat  plate Ss  we  have  from  Eq,  (ij..l) 
the  following  similitude  involving  tip-bluntness  and  boundary  layers  For  similar 
bodies  generated  from  the  same  equation 


the  following  correlation  applies s 


M 


KeA, 


/y  <  A> 


(Iul5) 


where  -f  represents  the  functional  dependence  concerned.  Two  alternative  forms 
of  Eq,  (I|,l5)  can  be  found  in  Appendix  I, 

It  is  essential  to  realize  that  the  above  correlations  do  not  imply  the 
possibility  of  obtaining  the  full  similitude  for  the  entire  flow  fields  which 
should  then  include  the  density  and  temperature  fields  within  the  entropy  layer. 
The  latter  flow  quantities  in  the  entropy  layers  have  been  excluded  from  the 


lil 


present  framework  of  study ^  however »  These  correlations  based  on  the  2.ero=order 
theory^  Eq.  will  be  subject  to  the  same  errors  as  Eq.  (I;,!)  j  ioeo 


0 


The  most  serious  among  these  is  presumably  that  associated  xd.th  </'  o  Implicit 
in  the  above  results  are  also  the  assumptions  of  the  linear  viscosity-temperature 
relation,  unit  Prandtl  number,  as  well  as  uniform  T„/t,  •  Thns  foi-  a  more  ac¬ 
curate  correlation  under  more  general  conditions,  the  similitude  based  on  our 
zero-order  theory  may  not  suffice.  As  a  part  of  the  study  complementing  the 
foregoing  analysis,  an  examination  of  the  hypersonic  similitude  under  less  re¬ 
strictive  conditions  will  be  given  in  Section  5<> 


SECTION  5 


THE  VISCOUS  HYPERSONIC  SIMILITUDE  INVOLVING 
LEADING=EDGE  BLUNT  NESS 


For  both  experimental  investigations  and  practical  applications j,  laws  of 
similitude  for  correlating  flow=fie.ld  properties  under  different  free-stream 
conditions  and  body  dimensions  are  desirable »  The  similitude  arrived  at  pre™ 
viously  in  Eq,,  (l+.l^)  is  simple  but  restrictive »  In  the  followingj,  the  more  gen¬ 
eral  viscosity-temperatui-e  relation  ^  assumed^  the  Prandtl  number  is 

taken  constant  but  otherwise  arbitrary j,  and  the  body  surface  temperature  is  arbi¬ 
trary,,  A  perfect  gas  with  constant  specific  heats  is  stall  employed  but  the 
specific  heat  ratio  is  no  longer  required  to  be  close  to  unity. 

The  inviscid  hypersonic  similitude  which  takes  into  account  the  effect  of 
slight  blunting  has  been  given  by  the  author^^.  In  a  recent  paper^,  Hayes  and 
Probstein'^  consider  the  similitude  of  hypersonic  boundary  layers.  From  these 
two  studies  a  viscous  hypersonic  similitude  involving  bluntness  is  not  difficult 
to  infer,  Hayes  and  Ihaobstein  cautiony  however^  that  the  similitude  involving 
bluntness  may  not  be  extended  to  the  viscous  case,,  or  only  under  a  condition  so 
restrictive  as  to  make  the  result  unlnterestir^.  In  view  of  thisy  an  examination 
of  the  basic  requirements  which  underlie  hj^personic  similitude  seems  proper  at 
this  point, 

Basle  Requirements  for  the  Validity  of  hypersonic  Boundary-Layer  Similitude 
Involving  Tip  Bluntness 

Hayes  and  Probstein  obsexve  that  the  hypersonic  small-disturbance  theory 
necessarily  falls  in  the  vicinity  of  the  blunt  nose  itself,,  and  it  would  appear 

A  certain  part  of  the  study  presented  in  this  Section  is  based  on  unpublished 
work  on  the  subject  which  was  sponsored  ly  the  U,  S,  Air  Force  Office  of 
Scientific  Research  (Mechanics  Division), 


h3 


that  the  theory,  hence  the  "equivalence  principle"  deduced  from  it,  would  break 
down  in  the  layer  of  fluid  of  increased  entropy  nearest  the  body«  It  may  then 
follow  that  hypersonic  similitude  cannot  hold  within  the  entropy  layer  as  well 
as  the  boundary  layer  beneath  it<,  The  major  conclusion  of  our  following  discus¬ 
sion  is  that  a  full  similitude  capable  of  correlating  the  density  and  temperature 
distributions  near  the  outer  edge  of  the  boundary  layer  is  not  possible ,  but 
correlations  for  the  pressure  and  velocity  fields,  as  well  as  the  density  and 
temperature  distribution  within  the  boundary  layers  proper,  do  exist. 

Let  us  first  examine  whether  similitude  can  exist  in  the  inviscid  flow  re¬ 
gion,  On  the  basis  of  our  discussion  in  the  earlier  part  of  Section  3,  we  see 
that  under  requirement  (A) ,  the  approximation  U.  '>6\J  holds  throughout  the  in¬ 
viscid  flow  region  including  the  entropy  layer.  The  applicability  of  the  inquire 
ment  (A) ,  and  the  degree  of  approximation  involved,  has  been  examined.  From  the 
fact  that  U.  u  we  can  conclude  that  the  pressure  ^  and  the  velocities  /lA' 
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and  U.  do  obey  the  equivalence  principle,  and  hence  the  hypersonic  similitude  . 
For  the  density  and  temperature  in  the  entropy  layer,  however,  the  "equivalence 
■’.principle"  may  not  be  strictly  applicable  since  these  quantities  may  depend  on 
the  entropy  determined  upstream  at  the  nose  shock. 

For  the  boundary  layer,  we  observe  that,  with  the  wall  conditions  specified 
the  solution  will  be  completely  determined,  by  the  three  quantities  at  its  outer 
edges  the  pressure  ,  the  total  enthalpy  H^(=  Ha>)  >  velocity 

iig  ,  To  be  sure,  for  determination  of  the  boundary- layer  solution,  the  speci¬ 
fic  knowledge  of  flow  density  or  temperature  at  the  outer  edge  is  actually  not 
necessaiy.  Following  essentially  the  same  procedure  as  Hayes  and  Probstein  in 
examining  the  system  of  equations  governing  the  boundary  layer,  we  find  that  the 


condition  required  for  the  existence  of  boundary-layer  similitude  is  that  the 
quantities  s  Hg.  and  U.^  must  obey  the  equivalence  principle  in  the  outer 

floWo  But^  under  requirement  (A) $  and  for  M  large,  this  condition  can  indeed 
be  satisfied.  In  particular,  we  have  ILg  U  >  A  similitude  can  therefore  ex¬ 
ist  in  the  hypersonic  boundary  layer.  ¥e  may  note  that  the  actual  magnitude  of 
the  error  in  the  heat-transfer  rate,  the  skin  friction,  and  the  displacement 
thickness  incurred  by  the  approximation  U  will  not  differ  a  great  deal 

O 


from  AU^/\J 


In  fact,  if  we  assume  local  similarity  (which  is  appropriate 


for  the  present  application)  the  corrections  to  these  quantities  will  be,  respec= 


tively,  -1/2,  -3/2  and  +l/2  of 


Au^/U 


¥e  may  hence  conclude  that,  for  a  perfect  gas  with  constant  specific  heats 
and  a  viscosity-temperature  relation  with  the  exception  of  the  dis¬ 

tributions  of  density  and  temperature  in  the  outer  flow,  the  classical  hypersonic 
similitude  of  Tsien  and  Hayes  can  be  extended  to  the  boundary  layer,  even  under 
the  influence  of  the  leading-edge  bluntness,  provided  the  flow  regions  imme¬ 
diately  downstream  of  the  leading  edge  are  excluded  and  the  flow  Mach  number  is 
sufficiently  high,  i.e. 

„/-/ 

xfe)  «l  M»< 


¥hen  the  bluntness  effect  is  not  present,  the  requirement  will  presumably  be 
less  stringent. 

3  + 

The  tip-bluntness  effect  introduces  an  additional  parameter  into 

the  boundary- layer  similitude  discussed  in  Ref.  55 •  The  similarity  laws  correla¬ 
ting  shock  shape,  surface  pressure,  surface  heat-transfer  rate,  and  skin  friction 
^in  plane  flows  may  be  written  in  a  form  comparable  with  Eq.  (i).o5)  for  the  zero- 
order  theory; 


¥e  see  that  the  siinplicity  gauned  by  the  zero-border  theory,  'which  assumes  a  very 
strong  shock  and  a  very  small  €  s  consists  essentially  in  the  elimination  of 
the  parameters  -y .  t/t  and  »  Alternative  forms  of  Eq,  (5.1)  are 

given  in  Appendix  II  as  are  the  similarity  laws  for  axisymmetric  flows. 

There  are  two  important  limiting  cases  of  the  similitude  given  above.  One 
is  the  case  of  negligible  boundary-layer  displacement  effect,  the  other  that  of 
a  very  strong  shock  wave.  For  the  first,  the  interaction  parameter 
will  become  merely  a  scaling  factor,  and  we  have 

A  discussion  of  the  second  case  follows. 

Hypersonic  Similitude  Invol'^ing  Strong  Shock  Waves 

For  this  case,  we  require  the  shocks  to  be  very  stror^  so  that  terms  of 

'2(0)  ^ 

may  be  neglected  in  comparison  with  unity,  BVom  the  vie'wpoint  of  practical  ap¬ 
plication,  consideration  of  this  special  fom  is  not  trivial.  The  only  essential 
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difference  in  basic  requirements  betiseen  this  similitude  and  that  of  the  zero- 
order  theory  is  the  condition  on  y  .  The  zero-order  theory  furnishes  results 
in  a  convenient  form  for  comparison  mth  experimental  data  correlated  by  this 
similitude* 

Under  the  assumption  of  a  strong  shock,  the  Rankine-Hugoniot  shock  condi¬ 
tions  may  be  simplified.  Using  a  suitable  affine  transformation,  we  arrive  at 
a  form  of  similitude  comparable  with  Eqs.  (Iiel5)  and  (^*1)? 

I —  —1  3  i 


(5c2) 


•with 


X.J'tL 


'w 


-‘(f) 


Gon^Jarir^  this  -with  the  more  general  sitiili-fcude  of  Eq.  (^*1),  "we  see  that  in 


or 


the  above  form  of  simili-fcude  the  dependence  on  the  parameter 
on  ,  becomes  very  weak  as  a  result  of  the  shock  being  very 

strong*  On  the  other  hand,  comparison  -with  the  zero-order  result,  Eq*  (ii«l5) , 
shows  that  as  a  consequence  of  relaxing  conditions  on  Pr  $  cO  j  and  6  ,  we 

have  now  to  include  as  invariants  the  parameters  R-  ,  cO  >  G  j  as  -well  as 
dependence  on  </  and  on  in  the  strong-shock  similitude 

cannot  be  too  strong,  however^,  since  both  parameters  are  associated  xd-th  the 


.  first-order  €  effect 


SECTION  6 


CONGEUSION 

A  theoretical  study  has  been  made  of  combined  leading-edge  bluntness  and 
boundary-lsyer  displacement  effect  in  hypersonic  flow  over  thin  bodies#  Through- 
out  the  investigation^  a  continuum  flow  of  perfect  gas  with  constant  specific 
heats  is  assumed.  An  approach  has  been  developed  based  on  a  flow  model  consist¬ 
ing  of  three  adjoining  regions;  an  inner  Isoninar  boundary  layer^  an  outer 
(detached)  shock  layer ^  and  between  these  two  regions ^  an  entropy  layer.  Solu¬ 
tions  for  this  flow  model  have  been  obtained  by  considering  the  mutual  interac¬ 
tion  of  these  three  regions  under  the  requirements  of  the  hypersonic  small- 
disturbance  theory  and  the  shock-layer  theory.  In  the  boundary  layer,  a  smaH 

I 

pertTorbation  procedure  in  6  =  — -y-  is  developed  liiich  gives  the  local  similarity 
result  as  a  leading  approximation.  In  the  entropy  layer,  the  low  density  leads 
to  the  simple  pressure-area  relation  Eq«  (3o7)«  In  the  detached  shock  layer, 
the  results  of  existing  shock-layer  theory  are  applicable*  The  zero-order  ap¬ 
proximation  in  this  approach  corresponds  to  the  classical  result  of  Busemann. 

In  the  theory  developed  for  the  leading  appro xiraation,  the  problem  of  hyper¬ 
sonic  plane  flows  over  a  thin  body  of  an  arbitrarily  specified  shape  is  reduced 
to  one  governed  ty.a  second-order  ordinary  differential  equation,  i.e.  Eq«  (l+.l). 
The  simplicity  of  this  equation  permits  solutions  to  numerous  problems  of  inter¬ 
est,  including  the  combined  effect  of  leading-edge  bluntness,  boundary-layer 
displacement,  as  well  as  surface  inclination.  Continuous  transitions  are  exhi¬ 
bited  between  limiting  solutions  -sdiich  agree  with  existing  blast  wave,  strong 
viscous  interaction,  and  hypersonic  wedge  theories  for  small  e  .  Application 
of  this  zero-order  theory  to  the  case  of  finite  ^  ,  namely  Y  -  I0I4.O,  shows  a 

rather  small  discrepancy  from  the  exact  solutions,  as  far  as  the  surface  heat- 
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transfer  rate  and  shock  shape  are  concernedo  One  interesting  feattire  fotmd  for 
inviscid  flow  over  a  blunted  wedge  or  cone  is  an  oscillatory  pressure  decay  over 
the  afterbody  surface#  The  resvilt  can  be  interpreted  as  one  limiting  form  of 
shock-layer  reimpingement  at  glancing  incidence#  It  may  also  be  related  to  the 
’’pressure  undershoot”  observed  in  experimental  studies  of  hypersonic  flow  over 
blunted  slender  cones,  and  could  be  meaningful  for  studies  of  bluntness  induced 
boundary-layer  instability# 

For  the  pressure  and  skin  friction,  a  (consistent)  first-order  e  correc¬ 
tion  to  the  zero-order  result  appears  desirable*  The  procedures  for  obtaining 
higher-order  approximations  have  been  formulated  in  Sections  2  and  3*  “While  ex¬ 
plicit  solution  to  the  boundary-layer  equations  for  the  hi^er-order  correction 
is  presently  not  available,  simple  and  effective  methods,  such  as  that  of 
Ref*  k9}  may  be  applied  for  this  puipose#  It  may  be  pointed  out  that  the  results 
obtained  herein  are  also  applicable  to  a  thin  airfoil  of  arbitrary  planform  via 
the  ’’strip  theory”#  The  consistency  of  this  strip  theory  with  the  present  ap¬ 
proach  has  been  pointed  out  in  Ref#  itl  and  discussed  more  fully  in  Ref#  56# 

Hypersonic  similitude  has  been  examined  in  Section  5  for  the  present  prob- 
Ifflii  ytieve  both  bluntness  and  boundary-layer  displacement  effects  are  present# 

The  classical  similitude  of  Tsien  and  H^es  can  be  extended  to  include  these 
effects  for  the  correlation  of  pressure  and  velocity  fields,  as  well  as  flow 
quantities  of  interest  within  the  boundary  layers  proper,  provided  the  vicinity 
of  the  small  blunt  nose  is  excluded#  For  flows  involving  strong  shock  waves,  a 
simplified  similitude  is  obtained#  The  zero-order  results,  Tihich  require  a 
st3X)ng  shock  and  also  small  C  ■>  may  be  readily  compared  with  experimental 
data  correlated  in  this  similitude#  For  consideration  of  experimental  verifica¬ 
tion  of  this  similitude  and  comparison  of  the  experimental  results  with  the  zero- 
order  theory,  the  reader  may  refer  to  Ref# 
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APPENDIX  I 


ON  THE  COMPLETE  SOLUTION  TO  Y(YY0"=  e  Jtt/z 


Throughout  the  study  presented  in  the  text,  it  is  tacitly  assumed  that 
the  solution  to 

~  X./'  fe  =  ^onsUnh  (1,1) 


/  iA 

behaves  like  ^  s  as  required  by  the  existing  blast-wave  theory 

This  implies,  for  the  special  case  '/ — >1  and  Y  =  0  ,  that  the  only  solution 

of 


Ye  (Ye 


(1.2) 


admissible  to  the  present  study  is  that  of  Eq<,  (Uo5)  of  the  text,  l«e» 


On  the  other  hand,  a  complete  solution  can  actually  be  obtained  for  Eq,  (l»2), 

which  not  only  includes  the  blast-wave  result  Eq.  (l.3)^  but  also  exhibits  the 

feature  of  a  "free  layer"  in  the  limit  of  — ^0  ,  It  thus  appears  that  by 

employing  the  complete  solution,  it  would  be  possible  to  describe,  in  analytic 
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terms,  the  transition  from  Lighthill's  free  layer  in  the  neighborhood  of  a 
blunt  nose  to  the  blast-wave  dominated  region  far  downstream.  However,  as  far  as 
application  to  the  present  problem  is  concerned,  it  is  also  essential  to  examine 
the  compatibility  of  this  solution  with  the  small-perturbation  requirement  (A) 
which  underlies  Eqs,  (l.l)  and  (1,2).  The  main  conclusion  of  the  following  dis¬ 
cussion  is  that,  unless  the  small- perturbation  requirement  is  removed,  or  relaxed, 
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the  improvement  gained  by  using  the  complete  solution  cannot  be  establisheds  and 
the  particular  solution  pertaining  to  the  existing  blast-wave  theory  appears  to 
be  the  only  valid  result  under  requirement  (A) -(C).  Below,  the  complete  solu¬ 
tion  to  Eq.  (lo2)  will  be  presented  first,  and  the  question  of  consistency  with 
the  small-perturbation  requirement  will  then  be  examined.  The  existence  of  an 
explicit  integral  form  of  Eq,  (1,2)  for  the  complete  solution  was  pointed  out  to 
the  writer  recently  by  Mr,  J,  P,  Guiraud  of  ONERA,  France, 

For  convenience,  we  shall  drop  the  subscript  e  from  ,  Since  the 
differential  equation,  Eq,  ( 1, 2),  does  not  contain  the  independent  variable  ^  , 

we  may  let  =  YY"*?  and  the  equation  can  be  integrated  to  yield 

=  eJ?>t  (Y+cJ 

that  is, 

YY'=/eM(Y*Q 

Integrating  once ,  after  separation  of  the  variable ,  one  has  the  complete  solu¬ 
tion 

For  large  i  Y  C,  ,  Eq,  (l,5>)  provides  the  blast-wave  result, 
ioe,,  Eq,  (1,3).  Now,  for  sufficiently  small  ,  Eq,  (1,5)  niay  be  developed 
for  <C<  C,  )  namely 


=  -f-  -/7Jt  (-ti  +  Cj). 
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¥e  see  that  the  constants  and  Q  have  to  be  related  as 

I  ^ 


Cz--i  cT//^ 


and,  what  seems  to  be  more  significant,  that,  for  small  ^  , 


Y  =  Mt  ■■■  ■ 


This  is  no  more  or  less  than  the  equation  of  a  free  layer.  Now,  if  the  domain 
of  validity  for  the  solution  in  the  blunt-nose  region  and  that  for  the  hypersonic 
small-distance  theory  should  join,  the  above  asymptotic  solution  for  small  ^ 
must  be  identified  as  Lighthill’s  free  layer  in  the  vicinity  of  the  leading  edge. 
Accordingly,  C,  belongs  to  the  order  (t/sA)  ,  so  that  the  equation  for  the 


free  layer  may  be  written  as 


CL  -  0(1) 


(I»6) 


The  above  relation  implies  that  the  shock  layer  detaches  from  the  body  at  =  0(t) 
with  0(l)  ,  The  complete  solution  may  now  be  expressed  as 

[y  +  (cvi/eA)  ^  -  3  J  ["Y  +  ^ 

(1,7) 

The  blast-wave  result,  Eq,  (1,3)3  would  be  obtained  from  Eq,  (1.7)  by  neglecting 
terms  of  the  order  (o^t'y/sAY)  , 

The  question  of  whether,  and  when,  the  more  complete  solution  Eq,  (I, 7)  is 
compatible  with  the  small-perturbation  requirement  (A)  will  be  examined.  In  this 
respect,  the  asymptotic  solution  for  large  )  ,  i.e,,  the  blast-wave 
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result,  has  been  studied  In  Section  2  of  the  text,  where  the  range  of  its  validity 
has  been  establishedo  ¥e  shall,  therefore,  be  concerned  with  two  remaining  ques¬ 
tions;  (l)  can  requirement  (A)  be  satisfied,  when  the  difference  between  Eqs. 
(1,7)  and  (l«3)  is  not  small,  and  (2)  when  the  blast-wave  solution  dominates,  does 
the  contribution  from  the  terms  belong  to  an  order  lower  than  that 

of  the  error  in  the  hypersonic  small-disttirbance  theory. 

The  first  question  amounts  to 

if  eAy/at  ^  0(0  3  ^ 

Since  -0(1)  >  y= 

=  °(f)  ^ 


(1^9) 


Hence , 


■f-l 


'cLY' 


0(el) 


Izd. 


=  0(0 


(I. 10) 


Therefore,  when  the  free-layer  character  begins  to  predominate  the  solution  for 
flat-plate  afterbodies,  the  small-disturbance  requirement  cannot  be  satisfied. 
The  second  question  concems 

if  €  yk-  i  ^  j 


2-tkI 

)  (d.  ' 

Z  \dY-  J 


KeJkY/  ' 

Since  sJkycot  »l  ,  Y  oj  (eJUtf^ 


(1,11) 


and 


diC 


Y 


\€AyJ 


-fdr 


(1,12) 


9) 


thus 


(1.13) 


Therefore j  the  terms  in  the  blast-wave  dominated  range  are  too  small 

to  be  considered  as  a  valid  improvement  of  the  small-disturbance  theory., 

Our  conclusions  ares  (1)  The  complete  solutionj  EqSo  (l.^)  or  (1.7) j  exhi¬ 
bits  an  apparently  realistic  behavior  of  the  detached  shock  layer  over  a  flat- 
plate  afterbody j  providing  a  transition  from  the  Lighthill's  free  layer  to  the 
blast-wave  solution!  (2)  Unless  the  small-perturbation  requirement  (A)  is  re¬ 
moved,  or  relaxed,  the  complete  solution  cannot,  however,  represent  an  improve¬ 
ment  over  the  existing  blast-wave  result j  (3)  The  solution  may,  perhaps,  find  a 
valid  application  to  the  corresponding  problem  in  the  "piston  theory"  where  the 
small-disturbance  assumption  is  not  required. 
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APPENDIX  II 


ALTERNATIVE  FORMS  OF  HTPERSONIC  SIMILITUDE  FOR  CORRELATING 
SHOCK  SHAPE,  DISPLACEMEWT  THICKNESS,  SKIN  FRICTION, 
SURFACE  PRESSURE  AND  HEAT-TRANSFER  RATE  IN  PLANE 
AND  AXISYMMETRIC  FLOWS 


On  the  basis  of  the  hypersonic  small=disturbance  theory  and  the  boundary- 
layer  theory  for  ideal  gases  with  constant  specific  heats,  constant  Prandtl's 
number  and  the  viscosity-temperature  relation  ^  oC  ,  similarity  can  exist 
among  flow  fields  around  thin  or  slender  bodies.  For  bodies  having  similar 
shapes  and  similar  surface  temperature  distributions,  i«e. 


one  may  thus  correlate  the  corresponding  shock  shapes,  surface  pressure  and 
heat-transfer  rate  as 

n\/\. 

where  and  ^  represents  arbitrary  functions  of  the  variables  involved,  dif¬ 
ferent  for  different  flow  quantities.  The  displacement  thickness  of  the  bound¬ 
ary  layer  ^  and  the  skin- friction  coefficient  can  be  correlated  in  the 

same  forms  as  for  the  shock  shape  and  the  coefficient  of  surface  heat- 

transfer  ,  respectively,  and  therefore  will  not  be  repeated  hereafter. 

It  may  be  useful,  for  application  in  experimental  studies,  to  consider  cer¬ 
tain  alternative  forms  of  the  similitude  cited,  A  few  of  these  forms,  with 
Eq,  (ll,l)  unchanged,  are  given  belows 
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(II. 2a) 


L 

jri 

(Mrr‘'‘  J 


)  P  )TI  ^  ^  ) 


hyVij 


5y-3j> 

1+30 


H 


r>TTf^ 


j 


(lIo2b) 


(II. 2c) 


(n,2d) 
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and 


where  ^  is  the  same  as  ^  defined  in  Section  ij.,  with  <6=  L  , 


The  forms  of  Eqs,  (II. 2a)  -  (ll,2d)  given  above  have  the  property  that,  for  cases 
involving  very  strong  shock  waves,  or  TI  will  drop  out  from  the  right  of 

Eqs.  (II, 2a)  -  (lIo2d),  and  the  total  number  of  the  similarity  parameters  can 
therefore  be  reduced  by  one.  The  forms  (ll.2b)  -  (lIo2d)  have  the  additional  ad¬ 
vantage  that,  for  strong  shocks,  small  e  ,  unit  Prandtl  number,  linear  repre¬ 
sentation  of  the  viscosity-temperature  relation,  and  uniform  ?  J^ot  only 

the  constants  Pr  and  <0  drop  out  from  the  similitude  but  also  the  parameters 
6  and  T^/%  are  not  required  to  remain  fixed  in  the  simulation.  The  last 
three  forms  correspond  to  the  forms  of  particular  solutions  presented  in  Figs.  3- 
6, 
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table  I 


COMPARISON  OF  SKIN  FRICTIONS  AND  SURFACE  HEAT-TRANSFER 
RATES  OF  A  FLAT  PLATE  AT  HIGH  MACH  NUMBER  PREDICTED 
BY  SOLUTIONS  ASSUMING  LINEAR  VISCOSITY-TEMPERATURE 
RELATIONS  AND  THOSE  BASED  ON  YOUNG ‘S  FORMULA 
FOR  ^ 


Table  1(a)  Asymptotic  Formulas  for  high  M 


Chapman  & 

Rubesin^^ 

Present 

Young^^ 

Reference  Temperature 

T„ 

1“  *  3  ^1 

^  1  To  J 

a-l 

©■ 

Ci-I 

1  1  Tw 

Z  2  To  J 

(P-)" 

r 

*  Wi\ 

2 

1 

1 

I'M 

m  x-T  Comparison  of  the  skin  friction  predictions  for  Pr  =  Oo71 

Table  1(b)  ^  0,76,  at  high  M  _ 


Jk. 

To 

(^f)  i'^f) 

Chapman  &  s  ; 

Rube sin  Present  Young 

0 

s  1„03  s  1 

l/io 

lolU  ;  lo06  s  1 

1/6 

lolO  s  1»05  t  1 

1 

0„99  s  1»05  s  1 

6U 


TABLE  II 


COMPARISON  OF  THE  ZERO-ORDER  SOLUTIONS  WITH  THE  SELF-SIMILAR 
SOLUTIONS  TO  THE  PROBLEM  OF  STRONG  SHOCK-BOUNDART  LAYER 
INTERACTION  OVER  A  FLAT  PLATE  FOR  /  =  l.i^O  AT  DIFFERENT 
WALL-TO-STAGNATION  TEMPERATURE  RATIOS 


Tw 

MYs 

Ay 

Sources  From 
Which  Data 

(0> 

'5 

H 

t; 

M 

5  3,, 

are  Given 
or  Deduced 

Y, 

■p 

“w 

Q| 

0 

0,62 

0,158 

0.0965 

0.216 

Whale n^^ 

0,82 

0,85 

0,89 

0,80 

(0.183) 

(0,10) 

(0.20) 

Lees^^ 

0,20 

-  - 

0,232 

0.118 

0,283 

Whalen^^ 

_  _ 

0,88 

0,90 

0.75 

O.liO 

_  - 

0,312 

—  - 

-  -  - 

Whale n^^ 

0.88 

-  _ 

1,00 

— 

0,5i;6 

_  »  _ 

0,855 

Stewartson^^ 

-  « 

0.89 

0.38 

NOTE; 

1,  Lees’  theory  based  on 

tangent-wedge  formula  for  pressure. 

2,  3  etc. 

represent  zero-order  solutions 

given  in  text,  .  , 

3«  Zero-order  theory  differs  from  others  only  in  assumption 

6  close 

to  one  (all  theories  assume  Pr 

®  1,  oC.  T  )  , 
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TABLE  III 


GOMPAEISON  OF  THE  ZERO-ORDER  SOLUTIONS  WITH  THE  SELF-SIMILAR 
SOLUTIONS  TO  THE  FLAT-PUTE  AFTERBODY  PROBLEM  FOR  /  =  l,iiO 


Ool5 

0,lj.0 

1.0 


V 

I  5 

/  aJ 

/O^ 

Sources  F3X)m 

r 

lo) 

£Xfrt 

7 

Which  Data  are 
Given  or  Deduced 

'<s 

■P^// 

0,73 

0.112 

.0671 

2.39 

Cheng  &  Pallone^ 

1.05 

i.Ui 

1.15 

(0.775) 

(0.121) 

(0.0697) 

(2.39) 

(Sakurai^^) 

(0.99) 

(1«33) 

(1.10) 

0.73 

0.112 

.0682 

2.56 

Cheng  &  Phllone 

1.05 

l.iUi 

l.lli 

n 

li 

.070U 

2.87 

It  It  ft 

It 

It 

1.10 

It 

.07W 

3.60 

It  It  It 

It 

It 

i.oU 

NOTES 


/Off  and  /o^  are  obtained  from  tte  self-similar  solutions  to  the 
boundaiy-layer  problem  under  the  pressiire  given. 


(for  y  “  O) 


The  numerical  coefficient  in  the  shock-shape  result  for  =  I.I4.O 

given  by  Cheng  &  Pallone  is  in  error;  the  corrected  formula  is 


-  0.73  . 


j  s  etc.  represent  zero-order  solutions  given  in  Eq.  (U.?) 

If  data  deduced  from  Sakurai’s  blast-wave  solution  are  used  as  a 
basis  for  comparison,  the  zero-order  approximation  will  generally 
appear  better. 

For  the  case  of  cylindrical  afterbodies,  the  blast-wave  solution  of 
Sakurai  gives 

Y,  =  0.795  -  0.0689  ,  and 

C/ Y,  -  0-958, 
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U  >>  a 


U  » 


ENTROPY  LAYER:  Vb  <  y  <  Ye 
SHOCK  LAYER:  Ye  <  y  <■  Ys 


(b) 


Figure  I  (a)  ILLUSTRATIONS  OF  THE  STREAM-LINE  PATTERN  OF  THE  FLOW  FIELD 
AROUND  A  SLENDER  BODY  WITH  SMALL  BLUNT  NOSE 
(b)  DIVISION  OF  THE  INVISCI D-FLOW  REGION  ABOUT  A  SLENDER 
AFTERBODY  INTO  A  SHOCK  LAYER  AND  AN  ENTROPY  LAYER 
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FIGURE  2  DETERMINATION  OF  THE  MAXIMUM  PERTURBATION  VELOCITY 
\lSu\^U^-cc  OVER  A  FLAT-PLATE  AFTERBODY  IN  AN  I NV I  SC  ID 
HYPERSONIC  FLOW 
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Ys/CKe’^)  _ 


10® 


10® 


(a)  SHOCK  SHAPE  (b)  SURFACE  PRESSURE 


10-^  I 

(c)  SURFACE  HEAT-TRANSFER  RATE 


Figure  3  RESULTS  BASED  ON  THE  ZERO-ORDER  THEORY  OF  A  FLAT  PLATE  UNDER  THE  COMBINED  EFFECT  OF  TIP  BLUNTNESS 
AND  BOUNDARY-LAYER  DISPLACEMENT  :{ REFER  TEXT  FOR  BASIC  ASSUMPTIONS,  RESTRICTIONS  AND  NOTATIONS) 


/  (?/ 


(a)  SHOCK  SHAPE 


(b)  SURFACE  PRESSURE 


(c)  SURFACE  HEAT-TRANSFER  RATE 


Figure  RESULTS  BASED  ON  THE  ZERO-ORDER  THEORY  OF  A  FLAT  PLATE  AT  INCIDENCE 
UNDER  THE  COMBINED  EFFECT  OF  BOUNDARY-LAYER  DISPLACEMENT  AND  ANGLE 
OF  ATTACK 


eAt 


Figure  5  RESULTS  BASED  ON  THE  ZERO-ORDER  THEORY  OF  A  FLAT  PLATE  UNDER  THE  COMBINED 
EFFECT  OF  LEADING-EDGE  BLUNTNESS  AND  ANGLE  OF  ATTACK 


Figure  6  RESULTS  BASED  ON  THE  ZERO-ORDER  THEORY  OF  A  BLUNTED  SLENDER  CONE 
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Attn:  Prof.  J.  Kaye  (1) 

Department  of  Aeronautical 
Engineering 
University  of  Michigan 
Ann  Arbor  j  Michigan 

Attn:  Prof.  Wo  C.  Nelson  (l) 

Therm  Advanced  Research 
Therm,  Incorporated 
Ithaca,  New  York 

Attn;  Dr»  A.  Ritter  (1) 

Engineeilng  Research  Institute 
University  of  Michigan 
Ann  Arbor,  Michigan 

Attn:  Profo  Oo  Laporte  (1) 

Prof.  Go  Uhleribeck  (1) 

Department  of  Astronomy 
University  of  Michigan 
Ann  Arbor,  Michigan 

Attn:  Profo  L«  Goldberg  (l) 

Department  of  Aeronautical 
Engineering 

University  of  Minnesota 
Minneapolis,  Minnesota 
Attn:  Profo  J.  D.  Akerraan, 

Chairman  (1) 

Institute  of  Meteorites 

University  of  New  Mexico 

Albuquerque,  New  Mexico 

Attn:  Profo  I.  LaPaz  (1) 


Institute  for  Mathematics  and 
Mechanics 

New  York  University 

25  Waver ly  Place 

New  York  3  j  New  York 

Attn;  Prof,  R.  Courant,  Dir.  (l) 

Guggenheim  School  of  Aeronautics 

New  York  University 

New  York  53,  Mew  York 

Attn;  Profo  J.  F»  ludloff  (1) 

Department  of  Mechanical  Engineering 

Northwestern  University 

Evanston,  Illinois 

Attn;  Profo  A,  B,  Cambel  (1) 

Department  of  Physics 
University  of  Oklahoma 
Norman,  Oklahoma 

Attn:  Prof,  R,  G.  Fowler  (1) 

Aerodynamics  Laboratory 
Polytechnic  Institute  of  Brooklyn 
Freeport,  L<.I«,  New  York 
Attn;  Prof.  A,  FerrL  (l) 

Palmer  Physics  Laboratory 
Princeton  University 
Princeton,  New  Jersqy 
Attn:  Prof,  W.  Bleakney  (1) 

Prof,  Wo  C,  Griffith  (1) 

Department  of  Aeronautical  Engineering 
Princeton  University 
Princeton,  New  Jersey 
Attn:  Profs.  S,  Bogdonoff, 

W,  Hayes  (1) 

Princeton  University  Observatory 

Princeton,  New  Jersey 

Attn:  Prof.  L.  Spitzer,  Jr.  (l) 

Department  of  Aeronautical  Engineering 
Rensselaer  Polytechnic  Institute 
Troy,  New  York 

Attn:  Prof,  R*  P,  Harrington  (1) 


EDUCATIONAL  INSTITUTIONS  (cont) 


Engineering  Center 

University  of  Southern  California 

University  Park 

Los  Angeles  7,  California 

Attnj  Dr»  R,  Chaun  (l) 

Guggenheim  Aeronautical  Laboratory 
Stanford  University 

Stanford,  California  (1) 

Defense  Research  Laboratory 
University  of  Texas 
P.  0,  Box  8029 
Austin,  Texas 

Attns  M,  J.  Thompson  (1) 

College  of  Engineering 
U*G«L«A« 

Los  Angeles,  California 

Attns  Dean  L.  M«  K»  Boelter  (l) 

Department  of  Physics 
U«G*L«A« 

Los  Angeles,  California 

Attn:  Poof.  J.  Kaplan  (1) 

Experimental  Research  Group 

University  of  Utah 

Salt  Lake  City,  Utah 

Attns  Profo  M.  A.  Cook,  Dir#  (l) 

Department  of  Chemistry 
University  of  'Wisconsin 
Madison,  'Wisconsin 

Attns  Prof.  J*  0.  Hirschfelder  (l) 

Sterling  Chemistry  Laboratory 

Tale  UnivQi’sity 

New  Haven,  Connecticut 

Attns  PTofo  Jo  Go  Kirkwood  (l) 


INDUSTRIAL  ORGANIZATIONS 

Aerojet  Engineering  Corporation 
6352  No  Irwindale  Avenue 
Box  296 

Azusa,  California  (l) 


ARO,  Inc. 

Tullahoma,  Tennessee 
Attn:  Dr.  R,  W.  Perry, 

Mr.  R.  Smelt  (l) 

AVCO  Manufacturing  Company 

lycoming  Division 

Stratford,  Connecticut 

Attn:  Dr.  J.  C.  Keck  (1) 

Convair 

San  Diego  Division 

San  Diego,  California 

Attns  Dr.  'W.  H.  Dorrance  (1) 

Fairchild  Engine  Division 
Fairchild  Engine  &  Aircraft  Company 
Farmingdale,  L.I.,  New  York 
Attn:  M:*s.  C.  Minck,  librarian  (l) 

General  Electric  Compary 
Research  Laboratory 
P.  0.  Box  1088 
Schenectady,  New  York 
Attns  Drs.  Nagamatsu,  White, 

Alpher  (1) 

The  Glenn  L.  Martin  Company 
Baltimore  3^  Maryland 
Department  ^20,  Mail  #3072 
Attns  Mr.  L.  Cooper  (1) 

Hu^es  Aircraft  Corporation 
Research  &  Development  Laboratory 
Culver  City,  California 
Attns  Dr.  A.  E.  Puckett  (1) 

Lockheed  Aircraft  Missile  Systems 
Division 

'y'an  Nuys,  California 

Attns  Dr,  L.  Ridenour  (1) 

Marquardt  Aircraft  Corporation 
7801  Havenhurst 
Van  Nays,  California 
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INDUSTRIAL  ORGANIZATIONS  (cont) 


Ifi-dwest  Research  Institute 

Department  of  Pl^sics 

I^0^^9  Pennsylvania  Avenue 

Kansas  City^,  Missouri 

Attn:  Mr*  K*  L»  Sandefur  (l) 

North  American  Aviation,  Inc* 
Aerophysics  Department 
122li|.  Lakemjod  Blvd* 

Domey,  California 

Attn:  Dr«  van  Driest  (l) 

Ramo-Wooldridge  Corporation 

8820  Bellanca  Avenue 

Los  Angeles  California 

Attn:  Dr,  M»  U»  Clauser  (l) 

Rand  Corporation 
1700  Main  Street 
Santa  Monica,  California 
Attn:  E,  P«  Williams, 

C*  Gaaley,  Jr*  (l) 

Sandia  Corporation 
Sandia  Base 

Albuquerque,  New  Mexico 
Attn:  Drs,  C.  C.  Hudson, 

M.  L.  Merritt, 

J*  D.  Shreve,  Jr*  (l) 

Stanford  Research  Institute 

Poulter  Laboratories 

Stanford,  California 

ilttn;  Drs.  Poulter,  Duvall,  Doll  (l) 


Mr,  J,  P.  Guiraud 
Irgenieur  de  Recherches  a’  1’ 
Oii.N«E*R*A» 

De  L* Office  National  D'Etudes  et  de 
Recherches  Aeronautiques 
29,  Avenue  de  La  Division  Leglend 
CHATILLON-sous-BAGNEUX  (Seine) 

France  (1) 


FOREIGN 


Institute  of  Aerophysics 
University  of  Toronto 
Toronto  Canada 
Attn;  Dr,  G,  N»  Patterson, 

Director  (1) 

Division  of  Mechanical  Engineering 
National  Research  Laboratories 
Ottawa,  Canada 

Attn;  Dr,  J,  Lukasiewicz  (l) 


